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PREFACE 

This  report  was  prepared  originally  as  a  series  of  lectures  given 
at  the  International  Summer  Course  on  '*Some  Aspects  of  Shallow  Water 
Oceanology”  held  at  Lunteren,  the  Netherlands.  The  subject  of  this 
phase  of  the  lectures,  "Generation  of  Waves  by  Wind,"  included  both 
deep  and  shallow  water  conditions.  The  decay  of  swell  in  both  deep  and 
shallow  water  was  also  discussed.  In  addition  to  the  original  prepared 
manuscript,  this  report  includes  some  material  resulting  from  the 
discussions  during  and  following  the  lectures. 
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GENERATION  OF  WAVES  BY  WIND 


STATE  OF  THE  ART 


I.  INTRODUCTION 


When  air  flows  over  a  water  surface  waves  are  formed.  This  is 
an  observable  phenomenon.  Just  why  waves  form  when  air  flows  over 
the  water  is  a  question  about  nature  which  has  not  yet  been  answered 
completely  or  satisfactorily  by  theoretical  means.  Furthermore,  why  do 
waves  have  the  heights  and  periods  that  are  observed?  All  theories 
begin  either  with  the  assumption  that  waves  do  form  when  wind  blows  over 
the  water  surface,  or  else  that  waves  must  already  exist  by  the  time  the 
wind  begins  to  blow,  or  else  the  theory  is  immobile. 

The  brothers  Ernst  Heinrich  and  Wilhelm  Weber  (1825)  were  the 
first  known  to  report  experiments  on  waves,  and  A.  Paris  (1871)  made 
actual  wave  observations  on  the  state  of  the  sea.  These  observations 
were  made  aboard  the  DUPLIEX  and  the  MINIERVA.  Although  Airy  (1848) 
did  theory  on  waves  and  tides,  wind  forces  were  not  included.  Other 
early  contributions  on  wave  observations  at  sea  included  Abercromby 
(1888),  Schott  ( 1893) ,  and  Gassenmayr  ( 1896) .  However,  the  best  early 
documentation  on  wave  observations  was  perhaps  that  prepared  by  Cornish 
(1904,  1910  and  1934).  Cornish  also  attempted  to  relate  wave  conditions 
to  meteorological  and  geographical  conditions. 

Rather  than  thinking  of  the  scientist  in  the  role  of  answering  the 
question  "Why  do  waves  form,  "  one  might  rather  think  of  him  as  a 
practical  engineer  who  knows  that  the  phenomenon  does  occur  and  who 
can  then  recommend  what  should  be  done  about  it.  Progress  is  made 
only  by  engineering  application  of  scientific  theory.  Theory  offers  no 
progress,  except  when  implemented;  otherwise  it  is  dormant. 


Since  a  little  theory  has  never  hurt  a  practical  engineer,  an 
oceanographer,  or  an  applied  scientist,  it  seems  quite  appropriate  to 
mention  various  theories  which  have  been  proposed  through  various 
stages  in  the  advancement  of  the  state  of  the  art. 


Although  Stevenson  (1864)  established  the  first  known  empirical 
formula  for  wave  generation,  the  classical  wotk  on  wind  wave  theory  was 
due  to  Lord  Kelvin  (1887)  and  Helmholtz  (1888).  The  Kelvin-Helmholtz 
theory,  which  can  be  found  in  Hydrodynamics,  by  Lamb  (1945),  pertains 
to  the  study  of  the  oscillations  set  up  at  the  interface  of  two  fluid  media 
of  different  densities  --  water  and  air,  for  example. 


Helmholtz  considered  the  media  of  mass  densities  and  P 

flowing  with  velocities  Uj  and  with  respect  to  each  other  as 
illustrated  in  Figure  1. 
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Figure  1 


The  interface  is  a  >vave 
and  velocity  of  the  upper  fluid.. 


‘urface,.  ^ 
and 


ar.d  I',  the  mass  densi'A 

A 

U2  the  mass  density  ^:rci 


velocity  of  the  lower  fluid,  respectively.  The  propagation 9i  veloc:*v  C 
is  the  speed  at  which  the  interface  travels  in  a  ic  rward  direction,  -.rd 
L  is  the  distance  between  two  successive  peaks  of  the  interface. 


Helmholtz  (1888)  showed  that  the  induced  o.scillatior,  if  smail 
compared  with  the  distance  L  took  the  form  of  a  wave  tram  st  *he 
interface  traveling  at  the  veiccity  C  such  that 

/’,(U,  -C)2.  ^  (.^)  (/>,■  Z’,' 


where  g  is  the  acceleration  cf  gravity  and  k  is  the  wave  numbe’ 
as  k  =  2  TT/L. 


Kelvin  (1887)  derived  the  same  resuJ’  m  a  d  iferent 
made  some  very  interesting  conclusions  For  example  wher  V 
it  can  be  shown  that 


C 
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A] 


If  one  considers  tne  upper  fluid  t-  be  -..ir  and  the  lowc  -  flu;d 
for  which  equal  to  about  1.  29  x  iO'  ‘  then  Eq  I  redi^ce - 

very  nearly  to  the  simple  form 


C 
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which  is  the  classical  equaticn  for  wave  celerity  obtained  from  ..r^e-r 
wave  theory. 

By  use  cf  the  quadra*ic  fv-rmula,  Eq  (1)  car  be  solved  :  a 

one  obtains: 


2 


c  = 


-rp-rr- 


A-/*. 

77^ 


A  A 


^2  -  «lj 


1/2 


(4> 


Now  one  can  write  Eq.  (4)  as  follows: 

c  =  TT+  tr 


(5) 


where 


and 


tr  = 


^2^2 

A  ^  . 


=  c 


yo  p 


^2>' 


(U^  -  Uj)' 


(6) 

(7) 


In  the  above  TT  and  “C  represent  an  average  of  the  corresponding 

2 

values  of  U  and  C,  and  is  the  expression  identical  to  Eq.  (2). 


If  is  less  than  the  term  involving  in  Eq.  (7),  it  will 

be  found  that  becomes  imaginary,  which  implies  a  condition  of  in¬ 
stability  in  the  development  of  the  waves,  and  this  leads  to  a  progressive 
increase  in  amplitude.  Under  these  conditions  the  wind  is  traveling  faster 
than  the  waves  and  there  will  be  a  continuous  transfer  of  energy  to  the 
waves,  which  in  turn  goes  into  the  form  of  increase  in  wave  height  and 
increase  in  wave  celerity.  The  term  -  Uj  represents  the  wind  velocity 
relative  to  the  water  and  is  usually  expressed  simply  by  U  =  • 

The  condition  of  instability  is  defined  when 


U> 


1.  //A 


(8) 


Since 


f !  f 

^  V  ^  z 


1.  29  X  10''\ 


one  may  obtain 


U  >  28  C 
^  o 


C 

o 


IT 


< 


”ZF 


(9) 


c 

The  ratio  has  been  defined  as  the  wave  age,  and  waves  are 

unstable  when  their  wave  age  is  less  than  1/28;  this  instability  manifests 
itself  as  a  progressive  increase  in  wave  amplitude. 
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Fmmn  <^#7  it  esm  ^  th&i  the  smallest 

velocity  that  a  capillary  wave  (such  as  a  ripple)  can  have  is  23.  2  cm/ sec., 
corresponding  to  a  wave  length  of  1.  7  cm.  and  a  period  of  .  073  seconds. 
Work  by  Crapper  (1957),  Schooley  (1960)  and  Pierson  (1961)  shows  that 
because  of  nonlinear  effects,  the  1.  7  cm.  is  somewhat  low.  According  to 
Kelvin  (1887)  the  waves  will  always  be  unstable  if  U  >  28  x  23.  2  cm/ sec. 
That  is,  if  U  >  6.  5  m/ sec  (12.  5  knots),  then  the  waves  are  unstable. 

U  =  6.  5  m/ sec  is  also  called  the  critical  wind  speed  required  for  gravity 
wave  generation.  According  to  Munk  (1947)  there  is  a  critical  wind  speed 
below  which  waves  do  not  form.  There  have  been  numerous  articles  on 
the  subject  of  critical  wind  speed,  some  supporting  and  others  objecting 
to  the  existence  of  a  critical  wind  speed.  Reference  is  made  to  the  work 
of  Cox  and  Munk  (1956).  Later  Munk  (1957)  appears  to  be  dubious  as  to 
whether  or  not  a  critical  wind  speed  exists,  citing  the  work  of  Mendel- 
baum  (1956)  and  Lawford  and  Veley  (1956).  If  a  critical  wind  speed  exists.- 
it  appears  from  all  literature  sources  that  it  exists  approximately  between 
2  and  6  meters  per  second. 

The  concept  of  critical  wind  speed  is  indeed  a  controversial  sub¬ 
ject  at  present.  Nevertheless,  there  is  still  belief  that  there  is  a  critical 
wind  speed  somewhere  between  4  and  6  m/sec.  ,  and  that  below  the  critical 
wind  speed  the  fluid  flow  is  hydrodynamic  ally  smooth  or  laminar  and 
above  the  critical  wind  speed  turbvdence  develops  and  the  interface 
becomes  hydrodynamic  ally  rough,  for  which  wave  amplitudes  increase 
with  time  and  distance. 

The  difficulty  with  this  theory  is  that  there  is  a  density  difference 
between  air  and  water  even  when  the  wind  does  not  blow,  and  yet  no  wavf  ^ 
are  formed.  Thus  density  difference  alone  is  insufficient  to  start  wa\-r 
generation;  the  waves  must  already  have  existed  by  some  other  means, 
then  they  can  propagate  as  free  gravity  waves.  If  there  is  a  critical  wind 
speed  the  wind  is  already  blowing,  but  there  are  no  waves.  With  an 
increase  in  wind  speeds,  waves  do  form.  Why? 

It  was  not  until  1925  that  Jeffreys  introduced  tne  theory  of  ' r  :r> 
hypothvasis,  "  based  on  the  concept  of  a  hydrodynamically  rough  re 
for  wave  generation.  In  his  paper  "On  the  Formation  of  Waves  by  Wind,' 
Jeffreys  (1925)  proposed  that  eddies  on  the  leeward  side  of  the  waves 
resulted  in  a  reduction  of  normal  pressure  as  compared  with  the  wind¬ 
ward  face  and  in  a  consequent  transfer  of  energv  from  wind  to  waves. 

His  results  suggested  that  the  wind  could  add  energy  to  waves  only  so 
long  as  the  wind  speed  was  equal  to  or  greater  than  the  wave  celerity, 
and  that  when  the  wave  celerity  became  equal  to  the  wind  speed  the  waves 
reached  maximum  height  and  the  sea  was  one  of  steady  state.  Also, 
lowest  wind  speed  required  for  wave  generation  was  on  ♦he  crde*'  o*'  two 
knots,  or  about  one  meter  per  second.  It  appears  that  the  value  :  tvA 
knots  is  the  best  accepted  value  for  the  critical  wind  speed. 

During  the  ten  years  following  the  work  of  Jeffreys,  somewhat 
more  detailed  observations  were  made  of  the  sea  surface  conditions. 
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For  example,  Schumacher  (1928)  reported  tlic  first  known  study  of 
"Stereophotography  of  Waves"  from  the  German  Atlantic  Expedition, 
and  Weinblum  and  Block  (1936)  also  reported  results  on  stereophoto- 
grammetric  wave  records.  This  latter  contribution  gave  results  of 
measurements  carried  out  on  board  the  motor  ship  SAN  FRANCISCO, 
under  observation  of  V.  Cornish.  Other  data  on  waves  during  this  period 
included  that  of  Williams  (1934),  who  reported  on  sea  and  swell  observa¬ 
tions,  including  early  methods  of  obtaining  data,  and  Whitemarsh  (1935) 
who  reviewed  data  on  unusual  sea  conditions  as  reported  by  mariners, 
and  discussed  the  cause  of  high  waves  at  -sea  and  the  effect  of  these  waves 
on  shipping. 

After  the  founding  of  the  Beach  Erosion  Board  in  the  War  Depart¬ 
ment  in  the  early  1930's,  serious  research  began  on  theory  and  forma¬ 
tion  of  gravity  waves.  The  first  important  report  was  completed  in  1941 
and  published  in  1948,  "A  Study  of  Progressive  Oscillatory  Waves  in 
Water,  "  by  Martin  A.  Mason  (1948).  This  report  updated  the  state  of  the 
art  to  about  1940. 

The  next  great  advance  in  the  theory  of  wave  generation  in  deep 
water  was  that  by  Sverdrup  and  Munk  (1947),  although  Suthons  (1945) 
had  already  prepared  forecasting  methods  for  sea  and  swell  waves. 
Whereas  .  '"freys  (1925)  took  into  account  only  the  transfer  of  energy  by 
normal  stress s^'s,  Sverdrup  and  Munk  considered  both  normal  and  tan¬ 
gential  stresses.  (See  Figure  2.) 

- 


Figv  re  2 

The  average  rate  at  which  energy  is  transferred  to  a  wave  by 
normal  pressure  is  equal  to 

=  -i:  (  p, 

where  w^  =  -  kAC  cos  k  (x  -  Ct)  is  the  vertical  component  of  the  particle 
velocity  at  the  surface,  and  p  is  the  normal  pressure  acting  on  the  sea 
surface.  L  is  the  wave  length. 
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The  average  rate  at  which  energy  is  transmitted  to  the  waves  by 
tangential  stress  is  equal  to 


1 

*17 


X 


a  Ox 
o 


(11) 


denotes  the  horizontal  component  of  particle  velocity  at  the  sea  surface, 
and  'X^  is  the  wind  stress. 


u  =  kAC  sin  k  (x  -  Ct)  (12) 

o 

^  ^  f  '  (13) 

2 

where  jf~  is  the  resistance  coefficient.  Various  experiments  and  obser¬ 
vations  have  been  made  leading  to  controversial  values  of  as  a 

function  of  wind  speed.  However,  a  number  of  different  authorities  appear 

to  have  advocated  a  value  of  jf  close  to  2.  6  x  10  ,  and  it  is  this  value 

utilized  by  Sverdrup  and  Munk  (1947). 

According  to  the  above  arguments,  the  energy  of  waves  can  in¬ 
crease  only  if  (Rj^  +  Rj)»  the  rate  at  which  energy  is  added  by  both 

normal  and  tangential  stresses  of  the  wind,  exceeds  Ry  ,  the  rate  at 

which  energy  is  dissipated  by  viscosity.  The  energy  added  by  the  wind 
goes  into  building  the  wave  height  and  increasing  the  wave  speed.  That 
is,  Rjj  +  R^  =  +  Rj^  ,  where  R^j  is  that  portiion  of  energy  trans¬ 

formed  into  wave  heights  and  R^  is  that  portion  of  energy  transformed 
into  wave  speed. 

During  the  early  stages  of  wave  develcpmert  mosr  of  the  energy 
is  transmitted  by  normal  stresses,  but  when  C/L  >  0.  37  the  trans¬ 
mission  by  tangential  stress  is  dominant.  The  effect  of  the  normal 
stresses  dominates  for  a  short  time  only.  During  the  time  that  the  waves 
are  growing,  the  effect  of  the  tangential  stress  is  most  important.  When 
C/U  =  1.0,  energy  is  added  jy  tangential  stress,  but  there  is  a  small 
amount  lost  due  to  normal  pressure,  and,  for  this  reason,  the  relation  is 
written  with  t  Rtvj  •  When  Ry  =  Rj  t  Rj^j  »  fbe  waves  are  said  to  have 

reached  maximum  height  and  celerity  fcr  a  particular  wind  speed  and 
are  independent  of  fetch  length  and  wind  duration.  This  condition  is  some¬ 
times  called  the  fully  developed  sea. 

According  to  the  work  of  Sverdrup  and  Munk  (1947),  the  solution 
of  the  hydrodynamic  equation  describing  wave  generation  entailed  a 
knowledge  of  certain  coefficients  or  constants  resulting  from  mathematical 
integration  which,  of  course,  could  not  be  determined  by  theory  alone. 

The  appropriate  constants  were  determined  by  use  of  empirical  data. 
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Hence,  there  was  no  way  of  knowing  whether  or  not  the  theory  was  correct. 
In  order  to  evaluate  these  constants  it  was  necessary  to  resort  to 
empirical  wind  and  wave  data,  which  at  that  time  was  very  limited.  Wind 
speeds,  fetch  lengths  and  wind  durations  were  estimated  from  meteor¬ 
ological  situations,  the  data  of  which  were  also  based  on  very  meager 
coverage.  The  waves  were  estimated  by  visible  means.  Out  of  this 
theoretical  investigation  grew  the  concept  of  the  significant  wave.  The 
significant  wave  height  was  estimated  as  the  average  wave  height  of  the 
waves  in  the  higher  group  of  waves,  which  later  became  identified  very 
closely  as  the  average  of  the  highest  one -third  of  the  waves  in  a  record 
of  about  20  minutes  duration.  The  significant  wave  period  was  the  cor¬ 
responding  average  period  of  these  waves. 

According  to  the  theory  as  evaluated  with  "ancient"  data  for  the 
significant  wave,  the  fully  developed  sea  resulted  in  the  following  relations; 

gH/U^  =  0.  26 


and 


=  C/U  =  1.37 

where  H  and  T  are  the  significant  wave  height  and  period  respectively, 
and  U  is  the  wind  speed.  It  then  became  quite  apparent  for  any  situation, 
either  wind  waves  or  swell,  that  a  whole  spectrum  of  waves  was  present, 
including  a  probability  distribution  of  wave  heights  and  a  probability  distri¬ 
bution  of  wave  periods.  Much  of  the  above  work  was  performed  during  the 
days  of  World  War  II.  Otherwise  earlier  publications  would  have  appeared 
in  the  literature.  In  fact,  as  early  as  1935,  the  Imperial  Japanese  Navy 
encountered  a  typhoon  in  the  Pacific  Ocean  and  many  observations  were 
taken  but  were  not  published  until  much  later  by  Arakawa  and  Suda  (1953). 

The  time  had  then  arrived  when  no  further  advance  in  wave  gen¬ 
eration  theory  could  be  made  without  reliable  recorded  data  and  an 
advance  in  statistical  theory  and  data  reduction  and  analysis.  There  is 
no  necessity  to  discuss  wave  recording  here  since  this  subject  is  well 
covered  by  Tucker  (1964). 

Barber  and  Ursell  (1948)  were  perhaps  the  next  to  present  a  very 
important  paper.  The  results  of  their  investigation  proved  the  existence 
of  a  spectrum  of  waves.  A  completely  new  field  of  theory  and  research 
had  been  initiated,  but  it  should  be  noted  that  oceanographers  were  slow 
to  take  advantage  of  this  concept.  This  research  had  laid  the  foundation 
upon  which  many  advances  have  been  made  in  the  "state  of  the  art, "  and  it 
is  because  of  this  research  that  the  Sverdrup-Munk  (1947)  works  are  con¬ 
sidered  "ancient." 

Although  much  research  was  carried  out  during  the  next  few  years, 
no  great  advances  in  the  state  of  the  art  were  published  until  Thijsse  and 
Schijf  (1949)  presented  wave  relationships  for  both  deep  and  shallow  water 
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ba^ed  on  wave  data  and  some  considerations  of  the  Sverdrup«Munk  theory. 
Eseperiments  on  a  paraffin  model  of  wind- generated  waves  by  Thijsse  and 
Schijf  show  a  high  negative  pressure  at  the  crest  of  the  wave*  which  is  in 
the  Sverdrup -Munk  concept  of  a  constant  wind  along  the  free 
surface*  but  which  is  in  accordance  with  Bernoulli's  equation  for  an 
increase  in  wind  speed  at  the  crest  axid  a  decrease  at  the  trough.  This 
experiment  was  certainly  a  great  contribution. 

Johnson  (1950)  applied  the  Pi -theorem  concept  for  dimensional 
analysis  and  presented  wave  relationships  for  deep  water  based  on  a 
collection  of  numerous  data  from  Abbots  Lagcon,  California.  At  the 
same  time  the  U.  S.  Army  Ccrps  of  Engineers  (1950)  had  presented  wave 
data  generated  under  hurricane  wind  conditions  icr  shallow  Lake  Okee¬ 
chobee,  Florida.  The  Corps  of  Engineers  also  presented  wind  and  wave 
data  for  inland  reservoirs.  Fort  Peck,  Montana  (1951),  and  later  for  Lake 
Texoma,  Texas  (1953).  Bretschneider  (1951)  presented  revised  wave 
forecasting  relationships  of  Sverdrup  and  Munk  ( 1947),  based  on  the  field 
data  of  Johnson  (1950),  laboratory  data  of  Bretschneider  and  Rice  (1951), 
and  numerous  other  data  collected  by  various  authors.  Bracelin  (1952) 
presented  an  unpublished  report  on  observing,  forecasting,  and  reporting 
ocean  waves  and  surf.  An  excellent  summary  of  wave  recordings  was 
presented  by  Wiegel  (1962). 

It  seems  that  the  year  1952  witnessed  the  first  acceleration  in 
the  "state  of  the  art.  "  Longuet -Higgins  (1952)  presented  the  Rayleigh 
distribution  for  wave  height  variability  based  upon  a  narrow  spectrum. 

Putz  (1952)  presented  a  Gamma-type  distribution  for  wave  height  and 
wave  period  variability  based  upon  analysis  of  25  twenty-minute  ocean 
wave  records.  Darbyshire  (1952)  and  Neum>tnr)  (  1952)  each  presented  wave 
spectra  concepts  and  relations  for  wave  generatior,  based  on  collection  of 
wave  data.  The  method  of  derivation  used  by  Neumar.’''  (1952)  is  contro¬ 
versial;  it  lead  to  the  introduction  cf  ?.  dimer sicnal  constart,  and  for 

high  frequency,  the  energy  was  found  to  be  pri.  portir '.a J  to  i  ,  where 

f  =  wave  frequency.  Watters  (1953)  derived  the  Hayleigh  distribution 

of  wave  height  variability  in  a  less  sophisticated  manner  than  Longuet- 
Higgins.  and  the  data  of  Darlington  (1954)  supported  the  Rayleigh  distri¬ 
bution.  In  fact,  the  Gamma -type  distribution  for  wave  heights  of  Putz 
(1952)  was  represented  very  closely  by  the  Rayleigh  distribution.  Ichiye 
(1953)  studied  the  effects  of  water  temperature  on  wave  generation. 
Homada,  Mitsuyasu  and  Hase  (195  3)  performed  laboratory  tests  rui  wind 
and  water. 

The  first  acceJeratior  of  the  "state  of  the  a-t"  did  net  rnea^  much 
in  regard  to  the  development  of  the  theory  of  wave  generation,  since  this 
was  purely  empirical,  including  both  I'ield  and  iiboratory  data  collection, 
except  for  the  work  of  Longuet -Higgins  (1952)  and  Watters  (1953). 

Ursell  (1956)  had  surveyed  the  problem  of  wind  wave  generation 
and  opened  with  the  statement  that  "wind  blowing  over  a  water  surface 
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generates  waves  in  the  water  by  a  physical  process  which  cannot  be 
regarded  as  known;"  he  concluded  that  "the  present  state  of  our  knowledge 
is  profoundly  unsatisfactory.  " 

Bretschneider  and  Reid  (1954)  presented  a  theoretical  develop¬ 
ment  for  the  "Change  in  Wave  Height  due  to  Bottom  Friction,  Percolation, 
and  Refraction;"  Bretschneider  (1954)  combined  these  relationships  with 
the  wave  generation  relationships  given  by  Sverdrup  and  Munk  (1947),  as 
revised  by  Bretschneider  (1951),  to  obtain  shallow  water  wave  generation 
relationships  for  wave  height  and  wave  period  as  a  function  of  wind  speed, 
fetch  length  and  water  depth.  Sibul  (1955)  investigated  in  the  laboratory 
the  generation  of  wind  waves  in  shallow  water.  Aside  from  the  above  refer¬ 
ences  and  the  contributions  of  Thijsse  and  Schijf  (1949)  and  the  U.  S.  Army 
Corps  of  Engineers,  Jacksonville  District  (1950),  there  is  still  another 
contribution  on  wind-generated  waves  in  shallow  water.  In  1953  and  1954 
Keulegan  performed  experiments  at  the  National  Bureau  of  Standarcs; 
but  as  far  as  is  known,  this  information  has  not  been  published.  However, 
it  has  been  ascertained  that  the  data  of  Keulegan  is  in  agreement  with  that 
obtained  for  Lake  Okeechobee  and  is  also  in  agreement  with  the  relation¬ 
ships  presented  by  Bretschneider  (1954).  Saville  (1954)  published  a 
useful  report  on  the  effect  of  fetch  width  on  wave  generation. 

Data  collection,  although  limited  in  quantity  and  quality,  also 
persisted  during  the  period  from  1950  to  about  1955.  For  example,  Unoki 
and  Nankano  (1955)  in  Japan  published  wave  data  for  Hachijo  Island; 

Titov  (1955)  presented  works  in  Russian;  and  Bretschneider  (1954)  did 
work  on  shallow  water  of  the  Gulf  of  Mexico. 

There  seems  to  be  a  transition  period  during  the  years  1955  to 
1960.  His  ideas  based  on  theory  and  verified  with  data,  Krylov  (1956  and 
1958)  presented  the  Rayleigh  distribution  for  wave  height  variability  and 
postulated  also  that  the  Rayleigh  distribution  applied  to  the  wave  length 
variability,  which  could  be  transformed  into  a  period  distribution  function. 
Bretschneider  (1957  and  1959)  also  verified  the  Rayleigh  distribution  for 
wave  height  and  wave  length  variability.  He  developed  a  distribution  func¬ 
tion  for  wave  period  variability  which  was  in  agreement  with  that  postulated 
by  Krylov  (1958)  and  was  in  very  close  agreement  with  the  Gamma-type 
distribution  function  for  wave  period  presented  by  Putz  (1952).  Roll  and 
Fischer  (1956)  made  a  revision  of  the  spectrum  by  Neumann  (1952),  elimin¬ 
ating  the  dimensional  constant,  and  found  that  the  energy  at  high  frequency 

-5  -6 

was  proportional  to  f  instead  of  f  according  to  the  spectrum  of 
Neumann  ( 1952) . 

It  then  became  apparent  that  for  a  narrow  spectrum  it  was  safe  to 
assume  that  the  Rayleigh  distribution  applied  equally  well  for  both  wave 
height  and  wave  length  variability,  the  latter  readily  transformed  into 
a  wave  period  distribution  function. 

Statistical  representation  of  the  sea  by  wave  spectra  concepts 
through  the  work  of  Tukey  (1959)  and  Blackman  and  Tukey  (1958)  was 
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greatly  promoted  by  the  staff  at  New  York  University,  particularly 
Pierson  and  Marks  (1952). 

Under  the  assumption  that  the  joint  probability  distribution  of  wave 
height  and  period  was  uncorrelated,  Bretschneider  (1957,  1958,  1959) 
proposed  a  development  of  a  wave  spectrum  concept.  There  seemed  to 
be  some  similarity  between  the  Bretschneider  spectrum  and  that  proposed 
by  Neumann  (  1952),  and  finally  the  form  of  Bretschneider's  spectrum 
resolved  as  the  proposed  spectrum  of  Pierson  (1964)  based  on  the  simi¬ 
larity  theory  of  Kitaigorodski  ( 1961). 

At  this  time  also  Bret  schneider  (1958^  again  revised  the  wave 
forecasting  relationships  for  both  deep  and  shallow  water.  The  practical 
graphs  for  wave  forecasting  are  given  in  the  revised  version  of  Beach 
Erosion  Board  Technical  Report  No.  4  (1961).  These  relationships 
presently  are  undergoing  a  further  revision  which  should  increase  the 
accuracy  of  wave  forecasts. 

However,  during  the  ten  years  preceding  about  1955,  most  of  the 
effort  was  devoted  to  analytical  expressions  and  little  to  theory  of  wave 
generation.  During  the  days  of  Jeffreys  (1925),  Sverdrup  and  Munk  (1947) 
among  others,  the  concept  of  wave  spectra  was  not  promoted.  However, 
because  of  the  wave  spectra  enlightenments  of  Barber  and  Ursell  (1948), 
Seiwell  ( 1948),  Neumann  ( 1952) ,  Darbyshire  (  1952) ,  Bretschneider 
(1957),  Burling  ( 1959),  and  Pierson  (1964),  among  others,  a  new  channel 
was  opened  for  wave  generation  theory. 

At  this  point  mention  should  be  made  of  the  great  contributions  on 
wave  theory,  wave  probability  distribution  functions,  and  wave  spectra 
proposed  by  Miche  (1954).  In  particular,  Miche  proposed  the  Rayleigh 
distribution  to  wave  steepness,  a  theory*  not  previously  proposed.  This 
distribution  function  should  have  «?  wide  application  for  engineering 
studies. 

It  was  not  until  Phillips  (1957)  and  Miles  (1957)  that  additional 
theoretical  concepts  were  developed.  Sverdrup  and  Munk  (1947)  con¬ 
sidered  that  the  wind  was  constant  in  velocity  in  crder  to  develop  their 
theory,  but  this  was  proven  otherwise  by  Thijsse  and  Schijf  (1949). 

However,  Phillips  (1957)  considered  the  fact  that  the  wind  was  rapidly 
fluctuating  about  some  mean  value.  It  is  very  true  that  winds  blowing 
over  water  do  not  consist  of  streams  of  air  in  steady  and  uniform  motion 
but,  rather,  of  an  irregular  series  cf  "puffs"  and  "lulls"  carrying  eddies 
and  swirls  distributed  in  a  disordered  manner.  The  atmospheric  eddies, 
or  random  velocity’’  fluctuations  in  the  air,  ate  associated  with  random 
stress  fluctuations  on  the  surlace,  both  pressures  (i.  e.  normal  stresses) 
and  tangential  stresses.  The  eddies  are  borne  forward  by  the  mean  velocity 
of  the  wind  and,  at  the  same  time,  they  develop,  interact,  and  decay,  so 
that  the  associated  stress  distribution  moves  across  the  surface  with  a 
certain  convection  veiccity  dependent  upon  the  velocity  of  the  wind  and 
also  evolves  in  time  as  it  mr  ves  along. 
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It  is  these  pressure  fluctuations  upon  the  water  surface  that  are 
responsible  for  the  early  generation  of  waves.  The  tangential  stress  is 
not  considered,  but  Phillips  (1957)  states  in  some  cases  that  the  shear 
stress  action  might  not  be  negligible.  The  theory  is  in  agreement  with 
wave  observations  during  the  early  stages  of  generation,  but  as  C/U 
approaches  unity  there  are  other  wave  generating  processes  to  take  into 
account,  such  as  sheltering  and  the  effects  of  variation  in  shear  stresses. 
Although  this  theory  tends  to  an  under -estimation  of  wave  heights  for 
C/U  close  to  unity,  it  may  be  considered  as  a  great  advance  in  wave 
generation  theory  insofar  as  the  initial  birth  and  growth  of  waves  are 
concerned.  A  very  important  aspect  results  from  Phillips  (1957)  baised 
on  dimensional  considerations;  i.e.,  ‘for  high  frequency  components' the 

energy  varies  as  f 

Miles'  theoretical  model  for  the  generation  of  water  waves  is 
based  on  the  instability  of  the  interface  between  the  air  flow  and  the 
water.  The  theory  of  Phillips  predicts  a  rate  of  growth  of  the  sea  pro¬ 
portional  to  time,  whereas  after  the  instability  mechanism  of  Miles  takes 
over,  the  rate  of  growth  becomes  exponential.  The  Phillips  model  is  an 
uncoupled  model  in  the  sense  that  excitation  (air  flow)  is  assumed  to  be 
independent  of  response  (sea  motion).  The  Miles  theory  represents  a 
coupled  model  in  which  the  coupling  can  lead  to  instability  and  consequent 
rapid  growth.  There  can  be  little  doubt  that  both  mechanisms  occur  in 
any  practical  situation.  At  some  frequencies  in  the  spectrum  the  uncoupled 
model  will  govern  and  at  others  the  instability  model  will  govern.  The 
work  of  Miles  (I960)  is  a  recognized  contribution  on  wave  generation 
theory. 


Ijima  (1957)  presented  an  excellent  paper  on  the  properties  of 
ocean  waves  for  the  Japanese  area  of  interest.  This  study  included  valuable 
information  on  wave  spectra  obtained  under  typhoon  conditions.  Also  a 
decided  difference  existed  between  wave  spectra  obtained  on  the  open 
Pacific  Coast  and  that  obtained  for  the  coast  of  the  Sea  of  Japan. 

Another  important  effort  for  obtaining  wave  spectra  was  that  con¬ 
ducted  by  members  of  the  New  York  University:  Chase,  Cote,  Marks, 

Mehr,  Pierson,  Ronne,  Stephenson,  Vetter  and  Walden  ( 1957) .  All 
embarked  upon  a  great  task  of  obtaining  the  first  directional  spectrum  of 
a  wind- generated  sea  by  stereophotographic  techniques,  although  Weinblum 
and  Block  (1936),  about  20  years  earlier,  carried  out  measurements 
entailing  stereophotogrammetric  reproduction  of  ocean  waves  on  board  the 
motor  ship  SAN  FRANCISCO. 

Sulkeikin  (1959)  presented  the  Russian  methods  of  forecasting 
wind  waves  over  water,  v/hich  followed  from  his  earlier  works  on  the 
theory  of  sea  waves  (1956).  Burling  (1959)  presented  a  spectrum  of  waves 
at  short  fetches  and  found  a  range  in  values  of  m  in  the  high  frequency 

f  of  the  wave  spectrum.  Korvin-Kroukovsky  (1961),  in  his  Theory  of 
Seakeeping,  summarizes  much  of  the  early  work  on  wind  wave  generation. 
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In  particular,  this  book  includes  a  very  comprehensive  bibliography  on 
waves  and  wave  theory. 

In  1961  an  International  Conference  on  Ocean  Wave  Spectra  was 
held  at  Easton,  Maryland,  the  proceedings  of  which  were  published  in 
1963.  *  This  conference  included  about  30  presentations,  plus  discussions, 
and  had  an  attendance  of  less  than  100  participants,  representing  a  very 
large  percentage  of  the  scien^^ists  and  engineers  in  the  world  who  have 
been  contributing  to  the  advancement  of  the  science  of  ocean  wave  spectra. 
It  can  be  said  that  this  conference  brought  the  "state  of  the  art"  up  to  date. 
Known  and  unknown  properties  of  the  frequency  spectium  of  a  wind¬ 
generated  sea,  by  Pierson  and  Neumann  (1961,  1963)  was  the  most  logical 
paper  to  lead  off  the  program.  Unless  one  understands  the  concept  of  a 
fully  developed  sea,  the  work  of  Walden  (1961,  196.3)  might  be  misinter¬ 
preted  since  his  data  were  for  very  short  effective  fetches.  However, 
this  difficulty  was  clarified  in  the  discussion  of  Bretschneider  (1961,  1963). 
Numerous  discussions  followed  and  the  program  continued  well  on  its  way 
throughout  the  four -day  period.  The  conference  produced  two  sources  of 
directional  spectrum:  Longuet-Higgins,  Cartwright  and  Smith  (1961,  1963) 
and  Munk  (1961,1963).  On  he  last  day,  with  heads  still  spinning,  it 
became  an  accepted  fact  that  the  one-dimensional  linear  concepts  were 
not  always  sufficient  to  describe  the  state  of  the  sea. 

However,  the  conference  was  not  intended  to  bring  forth  new  theory 
on  how  waves  form  when  wind  blows  ever  the  water,  except  for  discussion 
of  the  work  of  Phillips  (1957)  and  Miles  (1957),  and  an  introduction  "On 
the  Nonlinear  Energy  Transfer  in  a  Wave  Spectrum"  by  Hasselmann 
(1961,  1963).  Otherwise  the  theory  was  limited  to  that  required  for  data 
collection,  data  reduction  and  analysis,  and  dana  presentation  and  appli¬ 
cations.  The  work  of  Hasselmann  (1961,1963'  is  indeed  a  classical 
contribution,  but  it  still  doe  -  not  tell  us  why  waves  are  formed  according 
to  pre-described  elevations  and  irequercie?. 

During  the  final  stages  of  deveiopmer.t  of  a  wind -generated  sea, 
two  nonlinear  processes  could  become  signiizcart.  There  is  a  dissipation 
of  wave  energy  due  to  breaking  (whitecaps),  and  a  transfer  of  energy 
flux  between  frequency  bands  may  take  place.  The  theory  of  Hasselmann 
concerns  the  latter  process.  The  theory  shows  that  by  the  fifth  order 
interactions  energy  can  be  transferred  between  frequency  bands  in  the 
spectrum.  In  fact,  Phillips  (1960  and  1961,  1963)  has  shown  that  this 
theoretical  energy  transfer  can  occur  at  certain  third  order  resonant 
interactions.  So  far  spectral  energy  transfer  is  a  purely  theoretical  con¬ 
jecture  and  has  not  been  verified  by  observation  or  experiment. 


*  In  the  following  material  leferences  the  above  conferen*  e  and 
proceedings  are  shown  by  (1961,  1963), 
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Cartwright  (1961)  and  Pierson  (1961)  gave  brief  reports  on  the 
papers  presented  at  the  1961  conference,  prior  to  the  publication  of  the 
proceedings. 

After  the  conference  was  over,  the  participants  went  home  to 
work  again,  hoping  to  advance  the  state  of  the  art.  More  data  were  to 
be  collected,  and  this  required  further  development  of  instrumentation, 
and  an  advancement  of  statistical  theory  and  computation  procedures. 

After  the  conference  several  important  papers  appeared,  although 
they  had  probably  been  worked  on  fcr  several  years.  These  included  a 
paper  by  Korneva  (1961)  cn  wave  variability  which  tended  to  verify  the 
previously  mentioned  probability  distribution  functions;  a  joint  paper  on 
"Data  for  High  Wave  Conditions  Observed  by  the  OWS  'Weather  Reporter' 
in  December  1959"  by*  Bretschneider,  Crutcher,  Darbyshire,  Neumann, 
Pierson,  Walden  and  Wilson  (1962);  and  a  paper  by  Schellenberger  (1962) 
on  "under suchungen  uber  Windwellen  auf  linem  Binnensee.  "  Pierson  and 
Moskowitz  (1963)  contributed  a  new  form  of  the  one -dimensional  wave 
spectrum  based  on  the  similarity  theory  of  Kitaigorodski  (1961)  and 
found  out  that  this  spectrum  fell  somewhere  among  the  other  past 
proposed  spectra,  considering  the  inherent  errors  arising  from  diffi¬ 
culties  in  determining  and  defining  wind  speeds. 

Walden  and  Piest  (1961)  presente  data  and  analysis  on  wave 
spectra  obtained  near  the  Mellum  Plate  Ligi.ihouse,  located  ’n  the  some¬ 
what  sheltered  water  off  the  North  Sea  coast  of  Germany. 

A  good  summary  of  wave  theory  and  wave  generation  is  presented 
in  Volume  I  of  The  Sea,  edited  by  Hill  (1962),  particularly  Chapter  19, 
"Wind  Waves,  "  liy  barber  and  Tucker  ( 1962) . 

Numerous  data  reports  on  wave  spectra  are  now  becoming  available. 
For  exauiple,  Moskowitz,  Pierson  and  Mehr  (1962,  1963)  prepared  reports 
on  wave  spectra  estimated  froxn  wave  records  obtained  by  the  OWS 
"Weather  Reporter  I,  II  and  III"  and  the  OWS  "Weather  Explorer,"  and 
Pickett  (1962)  presented  wave  spectra  fcr  the  Argus  Island  tower  off 
Bermuda. 

Bretschneider  (1962)  presented  a  concept  on  modification  of  wave 
spectra  over  the  continental  shelf,  and  Ijima  (1962)  presented  an  interesting 
development  of  the  correlation  between  wave  heights  and  wave  periods  for 
shallow  water.  Kitaigorodski  and  Strekalov  (1962,  1963)  presented  con¬ 
tributions  to  an  analysis  of  the  spectra  of  wind  wave  motion  based  on 
experimental  data.  This  work  was  a  continuation  of  the  work  of  Kitai¬ 
gorodski  (1961). 

Goodknight  and  Russell  (1964)  presented  the  first  data  on  large 
wind  waves  in  shallow  water  of  the  Gulf  of  Mexico.  These  waves  were 
generated  under  hurricane  wind  conditions.  The  statistical  analysis  of 
the  data  showed  that,  for  all  practical  purposes,  the  Rayleigh  distribution 
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was  quite  satisfactory  for  representing  wave  height  variability  for  large 
hurricane  waves  in  shallow  water.  The  distribution  of  wave  periods  did 
not  follow  the  distribution  function  of  Bretschneider  (1957)  or  that  of 
Putz  (1952),  but  fell  somewhere  between  these  distribution  functions  and 
the  Rayleigh  distribution.  The  data  seem  to  consist  of  long  period  waves 
arriving  from  deep  water  combined  with  local  wind  waves  generated  at 
a  large  angle  to  the  swells. 

Hamada  (1963,  1964)  presented  twc  very  interesting  reports  based 
on  laboratory  experiments  of  wind  wave  generation.  The  f”^  law 
originally  j^iuposed  by  Phillips  il'*57)  was  stated  to  be  applicable  to  the 
limiting  boundary  of  the  instability.  For  very  short  fetches  and  high  wind 

-7  3 

speeds,  Hamada  (1964)  finds  the  high  frequency  relations  of  f  “  and 
f-8,94 


According  to  the  work  of  Bretschneider  (1959),  the  high  frequency 

m  2 

energy  varies  with  f  where  m  =  9  for  very  low  gF/U  ,  and 

2 

decreases  in  magnitude  to  m  =  5  for  very  large  gF/U  ,  corresponding 
to  fully  developed  seas.  Thus  there  is  agreement  at  initial  v^ave  genera¬ 
tion  for  f  ^  between  Hamada  (1963)  and  Bretschneider  (1959),  and  also 

_5 

agreement  at  fully  developed  wave  generation  f  among  Phillips  (1957), 
Bretschneider  ( 1959),  Pierson  ( 1963),  and  Hamada  ( 1964).  There  are 
still  efforts  required  for  the  exponential  part  of  the  spectral  equation, 

-  Bf 

i.  e.  e  where  B  is  a  constant.  According  to  Bretschneider  (1959) 

and  Pierson  (1963),  n  =  4  ,  but  the  factor  B  is  still  in  some  disagree¬ 
ment.  The  work  of  Burling  (1959)  indicates  that  n  might  be  larger  than 
n  =  4.  Bretschneider  (1961,  1963)  states  that  n  might  vary  between  4 
and  8  or  9. 

To  date,  all  theories  are  useful  in  attempting  to  understand  the 
mechanisms  involved  in  the  gene  ration  of  waves.  None  of  the  theories 
tells  us  why  waves  are  formed,  let  alone  why  the  wave  heights  and  periods 
are  as  observed  or  why  there  is  a  spectrum  of  waves.  However,  there 
is  enough  information  to  formulate  various  empirical  v/ave  forecasting 
relationships  for  certain  practical  applications.  The  accuracy  of  such 
wave  forecasting  relationships  depends  on  the  accuracy  of  the  wind  and 
wave  data  collected  and  used  for  the  empirical  relationships.  The 
accuracy  of  the  wave  forecasts  then  also  depends  upon  the  accurar\ 
of  the  wind  forecasts. 
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At  this  point  of  the  discussion  it  appears  in  order  to  introduce 
an  equation  which  describes  the  sea  state  wave  spectrunn,  including  the 
variability  of  wave  direction.  The  equation  can  be  written  as  follows: 


2  (g  A. 


mt  , 
e  -  1 

mt 


/oo 

77  (k.CJ  )  cos  (f> 


(14) 


y  U_cos  ^  ,  ilcji; 


g 


dT 


where 


m  = 


A  /  cos  jiS  '  ^ 

7^  i 


(J  ^ 


In  Eq.  (14)  TT  (k,  7!!! )  is  the  three-dimensional  pressure  spectrtim  as  a 
function  of  the  vector  wave  number  k  and  time  7^  ;  U  is  the  convection 
velocity  of  the  pressure  systems,  and  u,}c  is  the  friction  velocity  of  ±he 
shear  flow.  ^  is  the  coefficient  calculated  by  Miles  (1960),  and  / 
and  are  water  and  air  densities  respectively.  ^ 

There  is  little  need  to  extend  the  above  review  any  further  since 
the  directional  spectrum  has  been  discussed  quite  adequately  by  Tucker 
(1964).  It  is  hoped  that  most  of  the  important  contributions  on  wave  gen¬ 
eration  have  been  mentioned.  Eq.  (14)  represents  the  present  state  of 
the  art  on  wave  spectrum  generation  theory,  but  additional  empirical 
data  are  required. 


In  regard  to  practical  methods  for  wave  hindcasting,  Bretschneider 
(1964)  presented  a  paper  which  takes  into  account  the  complete  problem  of 
deep  water  waves,  storm  surge  and  waves  over  the  continental  shelf, 
the  breaking  wave  zone,  the  wave  run-up  on  the  beach  and  dunes  for  the 
March  5-8,  1962,  East  Coast  Storm.  This  paper  shows  the  results 
based  on  present  methods  of  wave  hindcasting  and  also  emphasizes  the 
areas  of  need  for  further  research.  A  very  important  consideration  of 
wind  wave  generation  over  the  shallow  water  of  the  continental  shelf  is 
that  of  the  total  water  depth.  The  total  water  depth  includes  the  combined 
effect  of  ordinary  tide  and  storm  surge.  The  various  problems  of  wind 
set-up  and  storm  surge  have  been  discussed  by  Bretschneider  (1958).  No 
further  discussion  on  tides  and  storm  surge  is  given  here  since  the  sub¬ 
jects  were  well  discussed  by  various  lecturers  at  Lunteren,  e.  g.  Drs. 

J.  R.  Rossiter,  W.  Hansen,  P.  Groen,  J.  Th.  Thijsse,  and  J.  C. 
Shonfeld.  Additional  work  on  storm  surge  problems  is  planned  for 
subsequent  reports. 
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Although  the  work  of  Sverdrup  and  Munk  ( 1947)  has  been  called 
"ancient,"  the  revised  wave  generation  parameters  are  still  considered 
the  most  practical,  and  the  application  of  these  relationships  still  results 
in  more  accurate  forecasts  than  the  other  underdeveloped  techniques. 

It  may  not  be  too  long  in  the  future  before  the  wave  spectra  techniques 
become  more  accurate. 
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II. 


PRACTICAi.  APPT  ^CATIONS  --  DEEP  WATER 


A.  SIGNIFICANT  WAVE  CONCEPT 

The  significant  wave  nnethod  of  wave  forecasting  was  that 
originally  introduced  by  Sverdrup  and  Munk  (1947)  and  is  sometimes 
considered  the  "ancient"  method.  The  wave  forecasting  parameters 
presented  by  them  evolved  fron^  theoretical  considerations,  but  the 
actual  relationships  required  certain  basic  data  fcr  the  determination  of 
various  constants  and  coefficient?.  Hence,  this  method  can  be  called 
semi -theoretical  or  semi-empirical. 

The  significant  wave  method  entails  certain  definitions.  The  sig¬ 
nificant  wave  height  is  the  mean  cr  average  wave  height  of  the  highest 
1/3  of  all  the  waves  present  in  a  given  wave  train.  The  significant  wave 
period  represents  the  mean  period  of  the  significant  wave  height.  It 
was  found  from  the  analysis  of  wave  records  that  the  significant  height 
is  nearly  equal  to  that  height  reported  from  visual  observations,  and  for 
this  reason  there  was  sometimes  a  certain  amount  of  agreement  between 
various  empirical  formulas  used  prior  to  the  development  of  the  theory. 

It  might  be  mentioned  that  the  significant  wave  period  represents 
a  period  around  which  is  concentrated  the  maximum  wave  energy.  From 
the  work  of  Putz  (1952),  Longuet-Higgins  (1952),  and  Bretschneider  (1959), 
the  distribution  of  the  various  wave  heights  can  be  determined  by  use  of 
the  significant  wave  height. 

The  wave  parameters  obtained  from  the  theoretical  work  of 
Sverdrup  and  Muiik  (1947)  car.  also  be  obtained  from  dimensional  con¬ 
siderations.  This  has  been  done  by  Johnson  (1950),  among  others, 
utilizing  the  Buckingham  Pi -theorem  (1914). 

The  factors  on  which  the  wind  wave  par' meters  for  deep  water 
depend  are  the  wind  velocity  U,  the  fetch  lengtn  F  ,  and  the  wind 
duration  t  .  Of  the  wave  parameters  only  wave  height  and  wave  period 
need  to  be  considered  since,  in  deep  water,  the  wave  length  L  = 

(g/2TT)T2  and  the  wave  celerity  C  =  (g/2T7)T.  The  waves  will  surely 
be  subject  to  the  influence  of  gravity  and  then  it  may  be  supposed  that: 

C  -  fj(U,F,t,g)  (15) 

and 


H  =  f2(U,  F,t,g) 


(16) 


17 


Eqs.  (15)  and  (16)  state  that  C  and  H  respectively  are  functions 
fj  and  of  U,  F,  t,  g,  but  apply  for  deep  water  only. 

From  the  symbols  appearing  in  Eqs.  (15)  and  (16)  one  may  write 
the  dimensions  for  deep  water  as  follows: 


Symbol 


Dimensions 


c 

LT 

H 

L 

U 

LT 

F 

L 

t 

T 

g 

LT 

Additional  relationships  can  be  written  if  the  water  depth  is  taken  into 
account. 


For  each  of  the  above  equations  there  are  five  variables  and  two 
dimensional  units,  whence  from  the  Buckingham  Pi-theorem  the  solutions 
will  each  be  functions  of  5  -  2  =  5  dimensionless  products,  with  2+1  =  3 
variables  to  each  product.  In  respect  to  Eq.  (15)  one  can  write 

O  =  Fj  (CU''*g^,  lU^'g^)  (17) 


With  respect  to  equation  (17)  one  <  b*3in=  the  dimensions 


and 


If  one  considers  wave  generation  in  shallow  water.  Eqs.  (15)  and  (16) 
become 


and 


C  (U,  F,t.d,  g) 

(15a) 

H  ::  (U,  F,  t,d-g) 

(16a) 

where  d  is  the  water  depth. 
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Equating  to  unity  the  sum  of  the  exponents  for  the  corresponding 
dimensions,  one  obtains  the  following  equations: 

1  +  a  +  b  =  0 
-1  -  a  -  2b  =  0 
1  +  c  +  d  =  0 
-c  “  2d  =  0 
e  +  f  =  0 
1  -  2  -  2f  =  0 


or 


The  simultaneous  solution  of  the  above  results  in 


a  =  -1 
b  =  0 
c  =  -2 


d  =  1 
e  =  -1 

f  =  1 


Using  values  of  the  above  exponents,  Eq.  ( 17)  becomes 


0 


C 

TT  = 


(18) 


In  a  similar  manner  the  corresponding  expression  from  Eq. 
becomes 


(16) 

(19) 


It  might  be  mentioned  that  the  Pi-theorem  is  a  most  powerful 
tool  if  properly  used*  It  is  extremely  important  to  realize  that  the 
expressions  for  physical  fact  must  be  dimensionally  homogeneous;  other¬ 
wise  there  are  some  scientific  factors  missing. 

Equations  (18)  and  (19)  represent  the  wave  generation  parameters 
for  deep  water,  based  on  dimensional  considerations.  \jj  ^  and  2 

are  functional  relations  that  must  be  determined  by  use  of  wave  data. 

gH/U  ,  C/U,  gF/U  and  gt/U  are  defined  respectively  as  the  wave 
height  parameter,  the  wave  speed  parameter,  the  fetch  parameter,  and 
the  wind  duration  parameter.  The  wave  speed  parameter  can  also  be 

written  »  which  is  a  better  form  because  the  wave  period  is  more 

easily  measured  than  the  wave  speed. 
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Using  the  above  parameters  Bretschneider  (1951)  revised  the 
original  forecasting  relations  of  Sverdrup  and  Munk  (1947),  utilizing 
much  additional  wave  data.  Before  1951,  however,  Arthur  (1947)  also 
revised  the  same  relations,  but  did  not  have  the  data  that  were  available 
in  1951.  Further  revisions  of  these  relationships  were  made  again  by 
Bretschneider  (1958).  These  forecasting  relationships  have  acquired 
the  name  S-M-B  method  for  Sverdrup,  Munk  and  Bretschneider. 


The  final  form  of  the  dimensionless  wave  generation  parameters 
appears  in  figure  3. 


The  curve  tU/F  is  the  relationship  between  wind  speed,  mini¬ 
mum  duration  and  fetch  length,  and  was  determined  by  numerical  inte¬ 
gration  of  the  following  relationships: 


(20) 


tU 

F 


Wave  forecasting  relationships  given  in  Figure  4  are  based  on  the  dimen¬ 
sionless  parameters  of  figure  3. 

For  long  narrow  bodies  of  water  such  as  man-made  reservoirs, 

rivers,  canals,  or  narrow  inlets,  corrections  need  to  be  made  for  the 

fetch  length.  Figure  5,  based  on  the  work  of  Saville  (1954),  can  be  used 

to  calculate  an  effective  fetch  length,  based  on  actual  fetch  length  and 

fetch  width.  The  effective  fetch  length  F  at  F  .  should  be  used 

°  e  min 

with  Figure  4  to  determine  significant  waves  for  long  narrow  bodies  of 
water. 


B.  COMPLEX  NATURE  OF  SEA  SURFACE 


The  significant  wave  description  is  a  simple  and  practical  means 
of  dealing  with  problems  in  wave  forecasting.  However,  it  is  important 
to  recognize  that  the  sea  is  very  complex,  made  up  of  many  variable 
heights  and  periods.  Figure  6  shows  a  schematic  interpretation  of  a 
typical  wave  record  which  might  be  obtained  from  a  wave  recorder.  The 
significant  wave  height  is  the  average  of  the  highest  1/3  of  the  waves  in  a 
given  wave  train,  of  at  least  100  consecutive  waves,  and  is  therefore  a 
statistical  parameter.  The  significant  wave  period  is  the  average  period 
of  the  highest  1/  3  of  the  wave  heights,  and  is  common  only  to  the  S-M-B 
forecasting  method.  The  Pierson-Neumann-James  (P-N-J)  forecasting 
method  considers  a  mean  apparent  wave  period  and  ranges  in  wave  period. 
Both  S-M-B  and  P-N-J  consider  the  probability  distribution  of  wave 
heights  and  in  both  methods  the  Rayleigh  distribution  is  used.  The  S-M-B 
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FIGURE  3  FETCH  GRAPH  FOR  DEEP  WATER 


FIGURE  4  DEEP  WATER  WAVE  FORECASTING  CURVES 

AS  A  FUNCTION  OF  WIND  SPEED, 
FETCH  LEN;''H,AND  WIND  DURATION 

(FOR  "-:'CHES  I  TO  1,000  MILES) 
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FIGURE  5  RELATION  OF  EFFECTIVE  FETCH  TO 

WIDTH— LENGTH  RATIO  FOR  RECTANGULAR  FETCHES 
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FIGURE  6  METHODS  OF  WAVE  RECORD  ANALYSIS 
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method  uses  a  probability  distribution  also  for  wave  period.  Wilson 
(1955)  introduced  the  space -time  concept  for  forecasting  waves  in  moving 
fetches  and  in  19('3  extended  the  work  for  use  on  a  high  speed  computer. 
Figure  7  is  reproduced  from  Wilr^on  (1955). 

Draper  and  Darbyshire  (1963)  presented  relatic  aships  for  fore¬ 
casting  maximum  wave  heights.  Figures  8,  9,  10  and  11  are  reproduced 
from  Draper  and  Darbyshire  (196  3)  where  it  should  he  kept  in  mind  that 
the  maximum  wave  is  equal  to  1.  6  times  the  significant  wave  height. 
Figures  8  and  10  are  for  deep  water  and  Figures  9  and  11  are  for 
coastal  waters. 

1.  Wave  Variability 


a.  Significant  Wave  Height.  The  significant  wave  height, 
as  mentioned  above,  is  a  term  common  to  both  the  S-M-B  and  P-N-J 
methods  of  wave  forecasting.  Just  how  the  significant  wave  height  is 
related  to  the  probability  distribution  and  also  the  wave  spectrum  can 
best  be  illustrated  by  figures  12  and  13.  Figure  12  shows  the  distribu¬ 
tion  of  wave  heights  as  visualized  by  a  histograph,  i.  e.  the  number  N 
(or  percent  P  of  the  total  number)  of  waves  in  each  wave  height  range. 
Figure  13  is  a  schematic  diagram  of  the  wave  spectrum  and  in  this  form 
is  called  the  period  spectrum. 


b.  Distribution  of  Wave  Heights.  The  significant  wave 
height  is  a  statistical  parameter,  and  about  16  or  17  percent  of  the  waves 
will  be  higher  than  the  significant  wave  height.  It  is  shown  by  Longuet- 
Higgins  (  1952),  Watters  (  1953),  Krylov  (1956  and  1958),  and  Darlington 
(1954)  and  verified  by  Bretschneider  (1959)  that  the  distribution  of  wave 
heights  for  a  narrow  spectrum  is  given  by  the  Rayleigh  aistribution.  The 
probability  density  for  wave  heights  is  given  by: 


p  (H)  dH 


-it,  e  ^  U|  dH 


(Zi) 


and  the  cumulative  distribution  is 

P  (H) 


1 


(22) 


where  H  is  the  individual  wave  height  and  H  is  the  average  wave  height. 


The  average  wave  height  H  ,  the  significant  wave  height  j 

and  the  average  of  the  highest  ten  percent  of  the  waves  Pjq  j  related 

as  follows: 
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FIGURE  8  GRAPH  RELATING  WAVE  HEIGHT 

TO  WIND  SPEED  AND  DURATION, 

AND  TO  FETCH;  FOR  COASTAL  WATERS 


FIGURE  9  GRAPH  RELATING  WAVE  HEIGHT 

TO  WIND  SPEED  AND  DURATION, 
AND  TO  FETCH;  FOR  OCEANIC  WATERS 


SURFACE  WIND  SPEED,  KNOTS  SURFACE  WIND  SPEED,  KNOTS 


FIGURE  10  GRAPH  RELATING 

WAVE  PERIOD  TO  WIND  SPEED 


AND  DURATION,  AND  TO 
FETCH;  FOR  OCEANIC  WATERS 


FETCH,  NAUTICAL  MILES 

(AFTER  DRAPER  AND  DARBYSHIRE  1963) 

FIGURE  II  GRAPH  RELATING 

WAVE  PERIOD  TO  WIND  SPEED 
AND  DURATION,  AND  TO 
FETCH  ;  FOR^^  COASTAL  WATERS 
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FIGURE  12  STATISTICAL 

DISTRIBUTION  OF  HEIGHTS 


FIGURE  13  PERIOD  SPECTRUM 
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The  cumulative  distribution  P(H)  is  given  in  Figure  14. 


c.  Distribution  of  Wave  Periods.  The  significant  wave 
period  is  not  a  statistical  parameter  of  the  wave  period  distribution  func 
tion.  The  significant  wave  period  has  statistical  significance  only  when 
both  wave  heights  and  wave  periods  are  considered.  Put?  (1952)  obtain.- 
a  Gamma-type  distribution  function  foi  wave  period  variability.  Krylov 
(1956  and  1958)  postulates  the  Rayleigh  distribution  also  for  the  wave 
lengths,  and  Bretschneider  (1959)  shows  that  the  distribution  of  the 
square  of  the  wave  periods  (proportional  to  the  deep  water  wave  lengths) 
can  be  represented  approximately  by  the  Rayleigh  distribution,  a  trans¬ 
formation  of  which  leads  to  the  following  distribution  function  for  wave 
periods; 


P  (T)  dT 


2.  7 


-0.  675 


e 


dT 


(24) 


and  the  cumulative  distribution  is 


-0. 675 

P  (T)  =  1  -  e 


(25) 


The  cumulative  distribution  P(T)  is  also  given  in  Figure  14.  Korneva 
(1961)  gives  results  which  verify  the  probability  distribution  functions  of 
wave  height  and  period. 

d.  The  Weibull  Distribution  Function.  In  many  cases  it 
may  be  desirable  to  represent  empirical  data  by  means  of  a  simple 
distribution  function.  Weibull  (1951)  proposes  a  simple  analytical 
distribution  function  for  use  in  certain  civil  engineering  problems  as 
follows: 


1  -  e 


-BX 


m 


(26) 


where  P  is  the  cumulative  distribution 
B  and  m  are  constants 

X  ^  0  is  the  variate,  i.  e.  H,  L,  or  T  as  the  case  may  be. 
The  probability  density  is  given  by 
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FIGURE  14  DISTRIBUTION  FUNCTIONS  FOR  PERIOD 

VARIABILITY  AND  HEIGHT  VARIABILITY 
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It  can  be  seen  that  when  m  =  2,  the  Weibull  distribution  function  is  the 
same  as  the  Rayleigh  distribution  (Eq.  21  or  22)  and  when  m  =  4  the 
Weibull  distribution  function  is  the  same  as  Eq.  24  or  25.  Thus  the 
Weibull  distribution  function  might  be  suitable  for  investigating  wave 
variability  when  the  data  deviate  somewhat  from  the  Rayleigh  distribution. 


where 


twice. 


The  moments  can  be  generated  from 


■Oo 


m 
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X^  p  (X)  dX  = 


1  \ii 

^  m 


(28) 


'o  \  I  \  I 

|~*  represents  the  Gamma  function. 

If  the  terms  of  Eq.  (26)  are  rearranged  and  the  logarithm  is  taken 
one  obtains 


In  Jn 


I  n  3  +  min  X 


(29) 


Eq.  (29)  is  that  of  a  straight  line.  When  In  In  is  .plotted  against 

InX  ,  the  intercept  becomes  in  B  and  the  slope  of  the  line  is  given  by 
m  .  P  is  the  cumulative  probability  distribution,  calculated  from  the 
data  and  B  is  related  to  the  mean  wave  height  or  mean  wave  period  as 
the  case  may  be. 


TT  "”2 

From  Eq.  (28),  m  =  0  ,  m,  =  X  ,  m^  =  X  ,  etc. 
follows  that  ° 
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It  then 


X  =  H  s  T  = 


1  \  na 

^1  r 


or 


B  = 
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1  +  -2- 
m 


1  +  -i- 

m 


m 


(30) 


Thus  Eqs.  (26)  and  (27)  become 


P  = 


(31) 
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P  (X) 


m 


m  - 1 


(32) 


By  use  of  Eq.  i  Kq.  (29)  becomes 


£.n  In  I  -^-'p 


m 


-£n 


m  ~ 


(33) 


It  must  be  emphasized  that  the  Weibull  distribution  should  be  used 
only  where  the  trend  of  the  dats  shows  a  nearly  line^.r  relationship  accord¬ 
ing  to  Eq.  (29)  or  Eq.  (3J).  If  such  a  linear  relationship  becomes  apparent 
from  the  data,  then  a  graphical  solution  is  possible,  and  a  mere  accurate 
solution  can  be  obtained  by  the  statistical  method  cf  least  squares. 

Although  the  Weibull  distribiition  has  no  theoretical  basis,  the  rancticn 
does  have  a  wide  range  for  practical  applications.  For  example.  Figures 
15  and  16  represent  ?.n  analysis  of  wave  data  according  to  the  Weibull 
distribution.  These  data  are  based  on  long  wave  records  obtained  at 
Lake  Texoma  by  U.  S.  Army  engineers  (1953)  and  are  summarized  in  a 
report  by  Bretschneider  (1959).  The  straight  lines  given  in  Figures  15 
and  16  were  fitted  visually,  although  a  least  squares  fit  might  have  been 
made  to  obtain  better  accuracy. 


e.  Joint  Distributioti  of  Wave  Height.s  and  Periods.  The 
joint  distribution  of  wave  heights  and  periods  is  quite  complex,  except 
for  the  special  case  of  zero  correlation  between  wave  height  and  wave 

period.  Figure  17  is  a  scatter  diagram  of  Tj  =  H/17  and  A  =  T  /  CT) 
for  a  case  very  nearly  zero  correlation.  If  the  number  of  wave  heights  is 
summed  independently  of  the  wave  periods,  one  obtains  the  marginal 
distribution  of  wave  heights,  and  if  the  number  of  wave  periods  squa’-^d 
is  summed  independently  of  the  wave  heights,  one  obtains  the  margir.a! 
distribution  of  wave  pc  riods  squared.  Beth  of  the  above  distributions 
are  represented  approximately  by  the  Rayleigh  distribution.  For  thr 
of  zero  correlation  between  wave  height  and  wave  period,  joint  di^t  - 
bution  is  given  directly  by 

p  (H;T)  -  p  (H)  •  p  (T)  ('--4) 

Based  on  Eqs.  (21),  (24),  and  (34),  the  joint  pimbability  c^stri- 
butions  have  been  calculated  for  variejus  ranges  of  H/H  and  T/T  and 
are  summarized  in  Table  I.  This  table  assumes  zero  correlation  between 
2 

M  and  T  Gumbel  (i960)  discusses  the  bivariate  exponential  distri¬ 
butions,  taking  into  account  correlation  coefficients  different  from  zero. 

Eq,  (34)  can  be  constdered  approximately  correct  for  pract:'  -1 
engineering  uses,  except  that  the  limit  of  breaking  waves  should  ai  >'s 
be  checked.  The  breaking  wave  limit  is  given  theoretically  by  Mic' 

(1954)  as  follows 


H 

L 


1 


tanh 


2TTd 

L 


34 


-2.0  - 1.0  0  I  2 

In  H 


FIGURE  15  SAMPLE  WEIBULL  DISTRIBUTION 
DETERMINATION  FOR  WAVE  HEIGHT 

(DATA  FROM  U.S.  ARMY  ENGINEERS,  1953) 
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2 


FIGURE  16  SAMPLE  WEIBULL  DISTRIBUTION 

DETERMINATION  FOR  WAVE  PERIOD 

(DATA  FROM  U.S.  ARMY  ENGINEERS,  1953) 


Relative  Wove  Height 


FIGURE  17  SCATTER  DIAGRAM  OF  »? 

AND  A  FOR  400  CONSECUTIVE 
WAVES  FROM  THE  GULF  OF  MEXICO 
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JOINT  DISTRIBUTION  OF  H  AND  T  FOR  ZERO  CORRELATION 


where  the  limit  in  shallow  water  is  given  by  the  solitary  wave  theory 


H 


0.  78 


(36) 


Based  in  part  on  theory  and  in  part  on  empirical  data,  Bretachneide r 
(1960)  presents  a  relationship  of  the  breaking  wave  limit  which  is  as 
follows: 


=  0.  124  tanh  — j- — 


1  +  0 


.152  tanh 


nd  ~ 


( 


Values  of  H  and  T  obtained  by  use  of  Table  I  should  always  bt 
checked  by  use  of  the  above  equations. 


C.  WAVE  SPECTRUM  CONCEPTS 


The  original  wave  spectrum  concept  for  forecasting  waves  is  that 
due  to  Neumann  (1952).  The  wave  spectrum  method  resulting  from  this 
concept  is  that  due  to  Pierson,  Neumann  and  James  (1955).  Much  work 
has  also  been  done  with  regard  to  a  wave  spectrum  method  by  Darbyshire 
(1952  and  1955).  Bretschneider  (1959)  has  added  a  wave  spectrum 
approach  to  the  significant  wave  method.  A  more  recent  spectrum  concept 
has  been  proposed  by  Pierson  (1964),  the  form  of  which  is  in  agreement 
with  that  derived  by  Bretschneider  (1959).  Whereas  the  significant  wave 
method  is  commonly  called  the  S-M-B  method,  the  Pierson,  Neumann 
and  James  method  is  commonly  called  the  P-N-J  method. 

Just  as  the  S-M-B  method  is  based  on  empirical  wave  data,  it 
is  with  the  P-N-J  method.  As  mentioned  before,  any  suitable  wave  1' re  - 
casting  method  must  be  calibrated  by  use  of  wave  data.  The  P-N-J  r.e^h  ri 
can  be  used  to  predict  the  spectrum  of  waves  from  which  one  may  ob^-i'n 
the  significant  wave  height  as  well  as  the  statistical  distribution  of  the 
waves.  The  S-M-B  method  is  used  to  predict  the  significant  wave  height, 
from  which  it  is  possible  to  obtain  the  wave  spectr  um  and  the  statistical 
distribution  of  the  waves.  Both  methods  utilize  the  distribution  function 
derived  theoretically  by  Longuet-Higgins  (1952).  This  distribution  function 
is  in  very  close  agreement  with  the  empirical  relationships  given  by 
Putz  (1952)  based  on  the  analysis  of  25  ocean  wave  records.  Consequently, 
when  both  S-M-B  and  P-N-J  methods  predict  exacUy  the  same  significant 
wave  height,  then  the  two  methods  result  in  exactly  the  same  distribution 
of  waves.  It  must  be  remembered  that  both  methods  are  based  on  actual 
wave  data.  However,  the  wu/e  data  are  not  the  same  utilized  ir.  .he 
two  methods,  and  the  method  of  analysis  of  the  data  is  different.  Tht 
analysis  of  the  data  for  the  S-M-B  method  consists  of  determinirig  the 
significant  wave  height  and  period,  which  in  turn  are  related  to  wind  speed, 
fetch  length,  and  wind  duration. 

Prior  to  the  development  of  the  P-N-J  method,  Neumann  (1952) 
proposed  a  theoretical  wave  spectrum  based  on  a  great  abundance  of 
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individual  wave  observations.  This  pectrum  should  be  considered  as  a 
semi -theoretical  or  semi -empirical  spectrum  in  the  strict  sense  of  the 
definition.  From  his  wind  and  wave  observations  Neumann  (1952)  computed 

the  parameters  H/T  and  T /V  ,  and  arrived  at  the  following  relation: 


const  e 


Irrv 


(38) 


(See  figure  18.  )^  The  parameter  H/'T^  is  dir  ec  tly  related  to  th^wave 
steepness  and  T/V  the  wave  age.  H  is  the  individaai  height.  T  is  the 
individual  wave  period,  which  he  called  the  apparent  wave  period,  defined 
as  the  time  between  two  successive  up-crc ssing-^  of  the  Stillwater  surface 
by  the  wave  form. 


From  this  empirical  relationship  (Eq.  38)  and  other  considerations, 
Neumann  (1952)  derived  his  theor  etical  wave  spectrum  of  energy.  From 
these  results  he  concluded  for  unlimited  fetch  length  and  duration  (fully 
developed  sea)  that  the  wave  height  was  equal  to  some  constant  times  the 

wind  velocity  to  the  5/2  power,  and  the  energy  at  high  frequency  was  f  ^  . 
Because  cf  the  method  of  solution,  a  dimensional  constant  was  introduced, 
but  the  equation  has  been  evaluated  by  considering  that  the  total  energy 
under  the  curve  of  the  spectrum  was  proportional  to  the  energy  of  the  equivalent 
root  mean  square  wave  height.  The  dimensional  constant  has  created 
objections  by  various  individuals  in  the  field  of  wave  forecasting,  particu¬ 
larly  those  who  adhere  to  dimensional  homogeneity.  However,  if  one 
accepts  the  relations  as  empirical  rather  than  theoretical,  then  these 
relations  ought  to  be  suitable  for  wave  forecasting,  at  least  over  the  range 
of  parameters  from  which  the  data  were  analyzed.  The  important  thing 
in  his  development  of  the  spectrum  is  the  fctct  that  tctal  energy  of  the 
spectrum  is  correct,  assuming  that  the  work  of  Longuet-Kiggins  is 
correct,  and  this  appears  to  be  the  case.  The  actual  distribution  of  the 
energy  with  respect  to  wave  pevi-d  or  frequency  might  not  be  exactly 
correct,  but  perhaps  is  in  close  agreement  with  the  actual  data  from 
which  it  was  derived.  Net  until  much  additional  data  become  available 
will  it  be  possible  tc  reconcile  any  differences  between  Neumann's  wave 
spectrum  and  a  mere  correct  spectrum. 


Roll  and  Fischer  (  1956)  suggested  a  different  derivation  of  the 
Neumann  spectrum.  In  this  derivation  there  resulted  no  dimensional 
constant,  the  wave  height  {(tv  fully  developed  sea  became  proportional 

2  -5 

to  U  ,  and  the  energy  at  high  frequency  was  f 

The  concept  of  the  wave  spectrum  is  certainly  a  great  advance  in 
trying  to  uiiderstand  the  nature  of  wave  generation.  It  is  shown,  for 
example,  in  spite  of  the  objections  to  the  Neumann  spectrum,  that  the 
waves  generate  from  the  high  frequency  end  of  the  spectrum,  and  with 
this  concept  Neumann  (1952)  derived  from  his  spectrum  relationships 
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FIGURE  18  RATIO  OF  WAVE  HEIGHTS  TO  THE 
SQUARE  OF  APPARENT  WAVE  PERIODS, 
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for  the  so-called  young  or  transient  state  of  sea.  That  is,  relations  called 
CO -cumulative  power  spectra  are  developed  from  which  it  is  possible  to 
predict  E-values,  where  E  is  related  to  the  generated  wave  energy. 

The  theoretical  wave  distribution  derived  by  Longuet-Higgins 
(1952)  is  directly  related  to  the  E -values,  as  this  condition  was  utilized 
in  the  development  of  the  wave  spectrum.  The  significant  wave  height, 
mean  wave  heignt,  etc.  are  related  to  E  as  follows: 


H 

ave 

1.  772 

X 

II 

2.  832 

(39) 

^1/10 

3.  600 

Typical  examples  of  wave  forecasting  relationships  based  upon  the  above 
concept  are  given  in  Figures  i9  and  20. 

The  spectrum  proposed  by  Bretschneider  (1959)  was  obtained  by 
use  of  a  theoretical  function  for  the  joint  probability  distribution  of  wave 
heights  and  periods.  The  period  spectrum  was  obtained  as  follows: 

S^(T)  =  f°H^  p  (H,  T)  dH  (40) 

■'O 

where  the  integration  is  over  all  wave  heights  as  a  function  of  wave  period. 
Assuming  no  correlation  between  H  and  T  : 

p(H,T)  =  p(H)  •  p(T)  ,  whence  (41) 


CO 


Sj,(T)  =  p(T) 


p(H)  dH  =  H^p(T)  =  ■^(ir)^p(T) 


By  use  of  Eq. 


(24)  one  obtains  the  period  spectrum; 


S_(T)  =  3.43(H)^— 2-t 
^  (T)^ 


675 

e 


For  the  frequency  spectrum: 

T  -  and  dT 


from  which  the  frequency  spectrums  becomes: 


SE(f) 


3.  43 


(H)^  ^-5  ^-:675  (fT)’^ 

0^ 


(42) 


(43) 


(44) 
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(  HO.  PUB  603,  1957) 


FIGURE  19  DURATION  GRAPH 

CO-CUMUL>t;vE  SPECTRA  FOR 
WIND  SPEEDS  FROM  20  TO  36  KNOTS 
AS  A  FUNCTION  OF  DURATION 


SISNIFICANT  IIIAVC  NCICMTUmI) 


(  H  0.  PUB  603 , 1957) 


RGURE  20  FETCH  GRAPH 

CO-CUMULATIVE  SPECTRA  FOR 
WIND  SPEEDS  FROM  20  TO  36  KNOTS 
AS  A  FUNCTION  OF  FETCH 
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SimiFtCANT  WAVE  NEWNTdMI) 


If  one  lets 


^  = 
U 


.625  gH33 


(45) 


gH  _  1  /  3 

TW  "  “rmr 


(46) 


it  then  follows  that 


Se(£)  =  3 


r  sr'i  1 

.43  - i-j 


«5 

f  ^  e  F 


.  675  27T  Uf 

IT  4  g 


(47) 


For  the  high  frequency  end  of  the  spectrum  Sj,(f)  is  proportional 

to  f”  .  This  is  in  agreement  with  the  results  of  Roll  and  Fischer  (1956) 
and  Phillips  (1957)  as  well  as  the  newly  proposed  spectrum  of  Pierson  and 
Moskowitz  (1963). 

Other  forms  of  the  wave  spectra  have  been  suggested  by  Bret- 
schneider  in  the  Proceedings  of  the  Conference  on  Ocean  Wave  Spectra 
(1963).  The  corresponding  period  and  frequency  spectra  are  as  follows: 

S(T)  =  a  t”^  e"^”^  (48) 


o/r\  ,-m-2  -bf' 

S(f)  =  a  f  e 


(49) 


For  the  Weibull  distribution  function  (m  =  n  -  1)  one  can  propose 
the  following  form  of  wave  spectra; 

,  -2liL  !(,/£  )-™ -1 1 

s(f/f)  =  ™  L  °  J  (50) 


where  the  peak  of  the  spectrum  occurs  at  f  =  f^  .  For  the  cumulative 
form  the  above  equation  becomes 


where 


S(f)  =  -  E  1^  -  e 
E  =  /' 

-'On 


(f/g""l 

m  o 

(51) 

S(f)  df 

(52) 

Pierson  and  Moskowitz  (1963),  based  on  the  similarity  theory  of 
Kitaigorodski  (1961)  propose  the  following  form  of  wave  spectra  for  a 
fully  developed  sea: 
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S(CJ)  dCJ 


(53) 


CJ 

where 

O)  =  2  rrf 

=  8.  1  X  10“^ 

^  =  .  74 

U,rt  c  '^a-s  the  measured  wind  speed  at  19.  5  meters  above  mean 
IV.  -> 

sea  level,  and  was  used  directly  instead  of  making  an  adjustment  to  the 
standard  anemometer  level  of  10  meters.  This  was  to  avoid  selecting 
someone  elses  drag  coefficient.  According  to  the  Weather  Bureau 
criteria,  the  wind  speed  over  water  at  elevation  19.  5  meters  is  about 
1.  1  times  that  at  elevation  10  meters. 

If  Eq.  (53)  is  integrated  to  obtain  the  total  energy  and  related  to 
the  significant  wave  height,  one  obtains 

gH/uf,  5  =  .21 

or 

gH/UjQ  =  0.254 
2 

which  is  very  nearly  equal  to  gH/U  =  .  26  used  by  Sverdrup  and  Munk 
(1947). 


Eqs.  (47)  and  (53)  are  directly  related,  any  differences  resulting 
only  from  the  empirical  relations  for  CJ ,  CJq>  ^2  functions 

of  wind  speed,  and  the  corresponding  interpretations  of  the  wind  speed. 

D.  FROUDE  SCALING  OF  THE  WAVE  SPECTRUM 

Assuming  that  the  Froude  law  applies  for  wave  spectra,  it  is 
possible  to  estimate  design  wave  spectrum  by  the  concept  of  Froude 
scaling.  The  design  wave  parameter,  for  example,  the  significant  wave 
height,  can  be  obtained  from  the  analysis  of  the  compiled  wave  statistics 
(wave  statistics  can  be  compiled  by  long  term  measurements  or  by  use 
of  wave  hindcasts).  Measured  wave  spectra  are  available  for  certain 
storm  conditions,  for  example,  Bretschneider,  et.  al.  (1962).  For  the 
design  storm,  once  in  100  years,  for  example,  the  significant  wave  height 
can  be  considerably  larger  than  that  corresponding  to  a  measured  wave 
spectrum. 
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In  a  report  by  Bretschneider  and  Collins  (1964),  the  most  seve 
hurricane  which  might  occur  in  the  Atlantic  Ocean  can  generate  a  62-i  .  t 
significant  wave  height.  The  measured  wave  spectrum  of  Bretschneidt- , 
et.  al.  (1962),  was  used  to  estimate  the  corresponding  wave  spectrum  f*r 
the  62-foot  significant  wave  height.  The  following  table  gives  signiiicant 
wave  heights  based  on  five  measured  wave  spectra.  I'or  conversion  to 
the  62-foot  significant  wave  height  the  Froude  length  scale  is  obtained 
from  A  -  62  (H  measured)  and  the  corresponding  time  scale  is 

'z  =JT. 


Table  II 


Date 

Time 

HS 

A  =  62  4  Hg 

17  Dec  1959 

0600 

35.  2 

1.  76 

1.  33 

1  f 

0900 

34.  4 

1.  80 

1.  34 

1  ! 

1500 

33.  9 

1.  83 

1.  35 

)  t 

1800 

39.  7 

1.  56 

1.  25 

18  Dec  1959 

0000 

35.  3 

1.  76 

1.  53 

To  convert  the  measured  frequency  spectrum  to  design  spectrum 

2  "s  6 ' 2 

the  ordinate  S(f)  ,  having  dimension  ft  sec,  .  must  be  multiplied  by  A 

and  the  abscissa  f  ,  having  dimension  sec  ,  must  be  multiplied  by  A 


To  convert  the  corresponding  period  spectrum  ^^o  de.sign  specfr 
the  ordinate  S(T),  having  dimension  ft  sec”  ,  must  ho  multiplied  bv 

X 


by  A 


and  the  abscissa  T  ,  having  dimension  sec,  must  be  multi pi:  -  a 
1/2 


A 


If  the  design  significant  wave  were  different  from 
and  X  would  change  accordingly. 


II  =  62  fees 


•  U  t-  ' . 


(Data  tabulated  by  frequency  analysis  is  often  given  for  intervals 
of  /Sk  f*  In  order  to  convert  such  data  to  a  design  .spect  rum,  one  must 

-1/2 

make  an  additional  conversion  by  multiplying  ^  f  A  .  ) 


For  a  number  of  reasons  the  period  spectrum  has  a  moic  pr.ic  ‘.u  i 
application  than  the  frequency  spectrum.  A  simple  operation  cor.vf  r*  ^  -h. 
frequency  spectrum  to  the  period  spectrum  as  follows; 


S(T)  dT  =  -  S(f)  df 


(64: 
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where 


f  =  1/T 
df  =  -  1/T^  dT 


Thus 

S(T)  =  ^  S(0  (55) 

Figure  21  shows  the  Froude  scaling  of  the  period  spectrum  for 

Hq  *  62  feet  as  given  in  Table  II.  The  Froude  scaled  spectra  in  Figure 
o  ( 

21  are  shown  superimposed  on  the  theoretical  spectri^m  given  by  Eq.  (43). 
The  significant  variability  at  some  frequencies  is  apparent.  The  theo¬ 
retical  spectrum  by  Eq.  (43)  can  be  considered  as  yielding  the  mean 
spectral  energy  density  in  any  frequency  band,  but  within  such  a  band  an 
individual  ordinate  will  vary  according  to  some  skew-type  distri¬ 
bution.  The  areas  under  all  spectra  given  in  Figure  21  are  identical  and 
are  related  to  the  same  value  of  the  significant  wave  height,  but  only  the 
individual  ordinates  vary.  For  example,  if  a  given  ordinate  varies 
according  to  the  Rayleigh  distribution,  then  the  upper  and  lower\95% 
confidence  levels  would  be  respectively  almost  twice  and  half  of  the 
mean  values. 

In  Figure  21  the  large  instabilities  in  the  scaled  observed  spectra 
at  low  periods  (high  fM^uencies)  are  not  quite  so  apparent  in  the  frequency 
spectrum.  These^nstabilities  are  probably  caused  by  the  introduction 
of  noise  in  the  original  record,  either  by  sample  rate  or  instrument  in¬ 
sensitivity,  or,  most  likely  a  combination  of  both.  Since  there  is  a 
smoothing  process  (related  to  A  f)  for  the  frequency  spectrum,  it 
becomes  apparent  that  a  similar  smoothing  process  (related  to  A  T) 
should  also  be  applied  to  the  period  spectrum.  This  is  because  A  T 
is  not  linearly  related  by  a  constant  to  1/Af*  Additional  studies  are 
required  to  investigate  the  disparity  observed  between  spectra  in  terms 
of  Af  and  those  in  terms  of  A 
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FIGURE  21  COMPARISON  OF  PERIOD  SPECTRUM 

FOR  62- FOOT  SIGNIFICANT 
WAVE  WITH  STATION  "j"  DATA 
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PROPAGATION  OF  WAVES  AXD  SWELLS  INTO 
SHALLOW  WATER 


When  waves  or  swells  ptcpagate  into  shallow  water,  a  number  of 
modifications  take  place:  refraction,  rhoaling  and  energy  losses.  These 
factors  will  be  discussed  subsequently ,  <  onsidering  only  a  simple  sinu> 
soidal  wave.  The  fact  th^.t  waves  tr  »vtl  in  groups  ol  variable  amplitude 
and  frequencies  --  the  wave  ■•pectrum  --  entail--  add;tional  considerations 
discussed  later. 

If  the  waves  are  long -ere -ted  and  n  :e  moving  obliquely  towards  a 
straight  shoreline  whose  depth  contours  are  also  straight  and  parallel  to 
^he  shore,  these  portions  of  the  wave  front  which  elfectively  feel  bottom 
first  are  retarded  first;  therefore,  the  wave  becomes  subject  to  a  pro¬ 
gressive  curving  or  refraction  which,  in  its  overall  effect,  tends  to  align 
the  wave  front  to  the  depth  contours.  This  is  in  analogy  to  simple  light 
rays  of  a  particular  wave  length  when  traveling  from  one  media  to  a 
media  of  different  density.  The  light  is  bent  or  refracted.  When  all  the 
wave  lengths  for  the  spectrum  ci  light  travel  from  cne  media  to  a  media  of 
greater  density,  each  wave  length  is  refracted  to  ^  different  degree,  which 
results  in  sorting  of  the  colors.  Similarly,  when  the  spectrum  of  ocean 
waves  enters  a  refracting  area,  the  different  frequencies  or  wave  lengths 
are  sorted,  and  the  resulting  spectrum  is  changed  accordingly.  Figure 
22  illustrates  the  effect  of  refraction  for  a  simple  wave. 

The  orthogonal s  (wave  rays)  represent  the  direction  that  the  wave 
fronts  are  taking.  They  become  curved  in  the  process  of  refraction  and 
in  general  may  tend  to  diverge  •;  r  converge.  Over  a  suamarine  canyon 
(Figure  23)  the  waves  will  always  tend  tc  diverge,  ■tr.d  conversely,  ever  a 
submarine  ridge  (Figure  24).  the  w^ves  wiil  always  tend  to  converge. 

It  is  generally  assumed  that  the  wave  energy  contained  oetween 
orthogonals  remains^constant  as  the  wave  front  progresses;  this  suppose^ 
that  there  is^  no^disper sion  r  f  energy  later-illy  along  ♦he  front,  no 
reflection  of  energy  from  the  rising  bettem,  and  none  lest  by  other 
processes.  If  b^  represents  the  distance  between  orthogonals  in  deep 

water,  and  b^  the  distance  between  crthogonals  somewhere  in  shallow 

water  (where  the  corresponding  wave  heights  are  and  ),  since 

the  energy  of  the  wave  is  proportional  tc  the  squ-.re  of  the  wave  height, 
it  follows  that 


o  c 


Cg, 


=  b  H  ‘ 

X  X 


or 


w. 


(56) 


(57) 
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FIGURE  24  CONVERGING 
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where  ,  the  refraction  coefficient.  The  above  is  a 

simplified  example  and  the  general  mechanics  of  transformation  for  a 
simple  wave  follow. 

Assuming  that  the  gradient  of  the  sea  bed  is  not  so  steep  that 
wave  energy  is  reflected,  and  that  there  is  no  gain  or  loss  of  energy  due 
to  lateral  difftic^ion  or  dispersion,  then  the  wave  energy  transmitted 
between  ortKo^nals  remains  the  same,  except  for  loss  due  to  bottom 
friction  and  percolation  in  the  permeable  bed.  Frictional  losses  represent 
the  energy  expended  in  overcoming  resistance,  shear  stresses  or  drag 
forces  on  the  sea  bottom.  Percolation  losses  represent  energy  lost 
through  actual  infiltration  of  the  wave  motion  into  the  semi-fluid  mass  of 
certain  finely  suspended  sediments  cf  the  permeable  sea  bottom. 

If  and  Pp  represent  the  frictional  and  percolation  energy 

losses  per  unit  area  of  bottom  per  unit  time,  respectively,  then  from 
Figure  25,  which  shows  the  elemental  length  ^  x  of  rising  sea  bottom 
in  the  direction  of  motion,  the  energy  loss  (time  rate),  in  terms  of  power 
(P)  transmitted  across  the  vertical  section  bh  between  the  orthogonals 
will  be: 

(P^  +  Pp)  b  ^x  =  [Pb  -  (5  (Pb)]  .  Pb  (58) 

whence 


-  =  -  (Pf  +  Pp)  b  (59) 

The  wave  energy  E  per  unit  surface  ar^^a  is  given  by 

E  =  /gH^  (60) 

( 

which  is  transmitted  across  the-  unit  area  with  tlie  grc  ip  velocity  C  ;  i.  e. 

g 

P  =  E  C  .  The  group  velocity  changes  from  deep  water  into  shallow 
g 

water  according  to 

C  =  nC  (61) 

where  C  is  the  wave  celerity  and  n  is  the  transmission  Coefficient 
and,  according  to  linear  wave  theory,  is  given  by 


i  [> 


2  kb 


where  k  =  2Jr/ir. — It  then  follows  that  the  power  of  transmission  is 


Pb  =  Eb  C  ^  g  b  (nC) 
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In  deep  water,  just  before  the  waves  are  due  to  feel  the  effects  of 
depth,  Eq.  (63)  becomes 


Pb  :.-lfgH^bnC 
oo  8y®o  co< 


The  ratio  of  Eqs.  (63)  and  (64)  results  in 


o  o 


Pb 

THT" 

o  o 


which  may  be  expressed  as 


^  K,p  (66) 

I 

where  K  ,  K  and  K-  are  known  as  the  refraction  coefficient,  the 
r  s  tp 

shoaling  coefficient,  and  the  friction-percolation  coefficient,  respectively^ 

I 

The  refraction  coefficient  is  separately  defined  as 


the  square  root  of  the  ratio  of  the  widths  between  orthogonal s  in  deep  and 
shallow  water.  According  to  linear  theory,  depends  only  on  wave 

period,  water  depth  contours  and  initial  deep  water  wave  direction.  In 
practical  coastal^nginecring  problems,  refraction  coefficients  can  be 
determined  by  one  of  two  methods:  wave  front  and  direct  ray  or  orthogcral 
methods,  both  of  which  are  discussed  in  B.  E,  B.  T.  R.  4(1961). 

The  shoaling  coefficient  is  given  compJcteiy,  according  to  linear 
wave  theory,  by:  _ 

/n”c  ' 


The  wave  celerity  C  in  shallow  water  is  related  to  in  deep  water  by 

C  =  tanh  kh  (69) 


and  in  deep  watfix— 1/2,  whence  from  Eqs.  (63),  (68),  and  (69 


K  =  1 

A 


Z  kh 


tanh  kh 


is  then  given  explicitly  as  a  function  cf  h  'L  or  h/L^  and  may  be 
obtained  from  tables  by  Wiegel  (1954)  where  =  H/H^. 
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The  friction-percolation  coefficient  K^p  is  defined  as 

/  V  1/2 

^  _  Pb  \ 

*^fp  -  'FT"" 

- -  0.0  / 

Returning  to  Eq.  (59)  and  substituting  Pb  from  Eq.  (71),  one  obtains 

d  (Kf/) 

P  b^  — —  =  -  (D,  +  D  )  b 
o  o  dx  £  p 


d  (K  ) 


(Dj  +  D  )  b 


According  to  Putnam  and  Johnson  (1949) 

^  _  4  „2  /fH^ 


°f  =  4  TT^ 


T-’  (sinh  kh)' 


and  according  to  Putnam  (1949) 


4  7T 


In  Eqs,  (74)  and  (75) 

=  mass  density  of  water 

f  s  dimensionless  friction  factor 

y  -  kinematic  viscosity  of  sea  water 

p  s  permeability  coefficient  of  Darcy's  Law  having 
the  dimensions  of  (length)^ 

Eq.  (75)  is  applicable  for  h/L  ^  0.  3  but  becomes  more  com* 
plicated  for  h/L  <C  0.  3. 

Returning  to  the  differential  Eq.  (73)  and  in  view  cf  Eqs.  '*4)  ano 
(75),  one  recognizes,  after  careful  examination,  that  this  is  of  the  •■xtrif 
general  form  as  Bernoulli's  nonlinear  differential  equation  rf  *he  Mrs* 
order:  i.  e.  I 


*  My  =  Ny" 
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Without  goirjg  into  the  aigebra,  one  may  v/rite  Eq.  (73)  as  follows: 


d 


F 


i 


and  the  solution  is: 


-  r  F  dx 
e 


(77) 


(78) 


In  general*  numerical  integration  cf  Eq.  (78)  is  required,  but  a 
number  of  special  cases  have  been  investigated  by  Bretschneider  and 
Reid  (1954).  Of  most  interest  at  present  are  two  simple  cases:  (1)  no 
bottom  friction,  and  (2)  no  percolation. 


The  first  case 


for  F,  =  0  results  in  K, 
f  fp 


K  =  e 
P 


-f 


F  dx 


v/here 


(79) 


The  second  case  for  F  -'0  results  in  K,  =  K,  where 

P  fp  f 


For  a  bottom  of  constant  depth  Eqs.  (79)  and  (80)  become 
respectively 


and 


X 


.  Fj  x] 


(80) 


(81) 


(82; 


Eqs.  (81)  and  (82)  can  be  used  for  practical  engineering  solutions  for  a 
bottom  of  variable  depth  where  small  increments  of  Ax  are  selected, 
each  increment  having  a  bottom  of  constant  depth  equal  to  the  mean  depth 
over  the  increment;  i.  e. 


where 


Kfp  ■ 

-F  Ax 

AKp  -  e  P  (84) 
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and 

[l.FjAx]"'  (85-, 

Graphs  presented  by  Bretschneider  and  Reid  (1954)  can  be  used 
to  facilitate  computations  by  the  numerical  method. 

Bottom  friction  and  shoaling  modifications  can  be  obtained  by  j-e 
of  Figures  26, and  27.  Figure  26  is  based  on  a  bottom  of  constant  depth, 
but  can  be  used  for  a  variable  bottom  by  representing  the  bottom  by  ^ 
series  of  sections,  each  having  an  average  bottom  depth. 


A  = 
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f  Ho  Kt  AX 
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FIGURE  26  RELATIONSHIP  FOR 

FRICTION  LOSS  OVER  A 

BOTTOM  OF  CONSTANT  DEPTH 
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FIGURE  27  K-  VERSUS  T /d 
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IV. 


GENERATION  OF  WIND  WAVES  IN  SHALLOW  WATER 


Less  information  avaiiible  on  wind  waves  in  shallow  water  than 
in  deep  water.  This  is  true  in  regard  to  both  theory  and  available  data. 

The  first  information  on  this  -ubject  was  given  by  Thijssc  {1949),  based 
on  limited  data.  Additional  date  and  relationships  were  brought  forth 
by  Dr.  Garbis  Keulegan  at  the  National  Bureau  of  Standards,  although 
never  published  to  the  knowledge  ol  the  author.  The  U.  S.  Army  Corps 
of  Engineers,  Jacksonville  District  (195S),  performed  an  extensive  field 
investigation  on  wind,  waves  at  '  tides  ^n  Lake  Oktec  icbee,  Florida. 

Based  on  the  hurricane  wind  and  wave  data  from  Li><t.  Okeechobee,  and 
some  ordinary  wind  wave  data  from  ‘he  shallow  rcf^/ior.r  ct  the  Gulf  of 
Mexico,  Bretschneider  (  1954>  was  able  to  esrablish  a  numerical  procedure 
for  computing  wind  waves  in  shallow  water  taking  bHom  friction  into 
account.  A  friction  factor  of  f  -  .01  appears  satisfactory.  Presently 
these  techniques  are  used  for  the  continental  shelf,  but  may  require  further 
calibration  when  more  wind  and  wave  data  are  available. 

A.  GENERATION  OF  WIND  WAVES  OVER  A  BOTTOM  OF 
CONSTANT  DEPTH 

If  d/T^  <C  2.  5  ft/sec^,  then  the  waves  effectively  "feel  bottom" 
and  the  depth  and  bottom  conditions  enter  as  additional  factors  with 
respect  to  the  heights  and  periods  of  waves  which  can  be  generated.  The 
effect  of  frictional  dissipation  of  energy  at  the  bottom  for  such  waves 
limits  the  rate  of  wave  generation  and  also  places  an  upper  limit  on  the 
wave  heights  which  can  be  generated  by  a  given  wind  speed  and  fetch 
length. 


Figure  4,  the  dimensionless  deep  water  wave  forecasting  relacior- 
ships,  in  effect  represents  the  generation  of  wave  e.iergy  in  deep  water 
as  a  function  of  F,  U  and  t  ,  since  the  energy  is  proportional  to 
whereas  Figure  26  represents  the  dissipation  of  wao  e  energy  due  to 
bottom  friction.  Figures  4  and  26  were  combined  by  a  numerical  method 
of  successive  approximation  tc>  obtain  relationship^:  for  the  generation  of 
waves  over  an  impermeable  bottom  of  constant  depth.  Best  agreement 
between  wave  data  and  the  numerical  method  was  cbt aired  when  a  bottom 
friction  factor  of  f  =  .01  was  selected.  Perhaps  a  "calibration  friction 
factor"  is  a  more  appropriate  term  since  it  would  t-'-.kt  into  account  other 
influential  factors  not  normally  i.ncluded  in  the  friction  factor  term  such 
as  energy  loss  by  "whitecaps.  "  The  curves  of  figure  28  are  the  results 

I  2  ^ 

of  these  computations.  The  curve  of  gT/U  versus  gd/U  is  based  or. 


2  2  2 

the  wave  data,  whereas  the  curves  of  gH/ U  versu-  gd/U  and  tzF/U 
are  based  on  the  numerical  computations.  The  curves  of  these  figures 
are  not  toe  much  different  from  those  presented  by  Thijsse  and  Schiji 
(1949),  reproduced  in  Figure  30.  Figure  29,  based  on  Figure  28,  gives 
wave  forecasting  curves  for  shallow  water  of  constant  depth  and  unlimited 
wind  duration  and  fetch  length.  Figure  28  may  be  used  when  both  the 
fetch  length  as  well  as  the  depth  are  restricted. 
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WIND  SPEED  IN  MILES  PER  HOUR 
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FIGURE  29  WAVE  FORECASTING  RELATIONSHIPS  FOR 
SHALLOW  WATER  OF  CONSTANT  DEPTH 
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growth  of  waves  in 


The  important  fact  from  the  above  material,  however,  is  the 
establishment  of  a  numerical  procedure  for  computing  wind  waves  in 
shallow  water  of  constant  depth  which  can  be  verified  by  use  of  wave 
data.  This  procedure  can  be  extended  to  a  bottom  of  constant  slope 
wherein  the  bottom  is  segmented  into  elements,  each  element  having 
a  mean  depth  assumed  to  be  constant.  Sample  computations  for  a  typical 
continental  shelf  are  given  by  Bretschneider  (  1957). 

Presently,  the  system  of  calculations  is  being  practiced  on  gen¬ 
eration  of  wind  wave  spectra  over  the  continental  shelf.  It  will  be  some 
time  in  the  future  before  the  system  is  perfected.  However,  Figure  31, 
reproduced  from  Bretschneider  (1964),  is  a  typical  example  resulting 
from  calculation  of  wave  spectra  in  shallow  water,  taking  into  account  the 
combined  effect  of  bottom  friction  and  wind  action.  The  hindcast  deep 
water  wave  spectra  is  based  on  the  method  of  Bretschneider  (1957)  and 
bottom  friction  is  taken  into  account  according  to  the  method  of 
Bretschneider  ( 1963)  wherein  the  wind  effect  is  added. 
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V. 


DECAY  OF  WAVES  IN  DEEP  WATER 


When  waves  leave  a  generating  area  and  travel  through  a  calm 
or  an  area  of  light  winds,  a  transformation  takes  place.  At  any  particular 
time  and  with  respect  to  decay  distance,  the  significant  wave  height 
decreases  and  the  modal  period  shifts  to  the  longer  period  waves, 
resulting  in  an  apparent  increase  in  significant  wave  period.  However, 
at  a  particular  decay  distance  the  modal  period  decreases.  To  under¬ 
stand  the  significance  of  a  period  shift,  one  must  refer  to  the  spectrum 
of  wave:  and  the  original  work  of  Barber  and  Ur  sell  (1948).  Recent 
work  of  Munk  (1964)  shows  that  most  of  this  transformation  takes  place 
in  the  first  1000  miles  of  decay.  Based  on  field  tests  of  waves  arriving 
from  the  southern  hemisphere  and  traversing  a  great  circle  to  Alaska, 
Munk  (1964)  shows  that  there  is  a  very  rapid  loss  or  transformation  of 
energy  during  the  first  1000  miles  and  a  shift  of  peak  frequency  to  about 
0.  6  to  0.  7  sec"^,  after  which  little  change  occurs. 

The  work  of  Wiegel  and  Kimberly  (1950)  and  Bretschneider  (1950) 
was  included  in  the  revisions  in  forecasting  decay  of  swell.  Empirical 
relationships  were  developed,  based  on  the  correlation  of  wave  data 
obtained  off  Southern  California  and  wind  and  fetch  characteristics  ob¬ 
tained  from  weather  charts  of  the  Southern  Hemisphere,  from  where  the 
waves  had  originated.  The  rather  meager  data  were  recorded  in  shallow 
water  and  had  to  be  transformed  for  deep  water  according  to  the  shoaling 
coefficients  given  by  Figure  27.  In  spite  of  the  above  difficulties, 
empirical  forecasting  curves  for  wave  decay  w'ere  determined  by  Bret¬ 
schneider  (1951)  and  are  presented  in  Figure  32.  The  scatter  of  data 
on  which  these  curves  are  based  indicates  that  the  above  relations  can 
deviate  by  as  much  as  t  50  percent.  Although  these  curves  are  not  com¬ 
pletely  satisfactory,  they  are  the  best  available  to  date.  Further 
refinements  are  certainly  necessary. 

Assuming  average  conditions  for  a  severe  storm,  a  wave  -pectra 
hindcast  was  made  for  the  end  of  the  fetch,  and  assuming  that  Figure  32 
applied  for  Hj^/H^  as  a  function  of  T  for  decay  of  the  wave  spectra, 

computations  were  made  for  various  elements  of  the  period  spectrum. 

The  results  of  these  calculations  are  shown  in  Figure  33.  The  curve  for 
D  =  0  is  zero  decay  distance  or  the  end  of  the  fetch.  Other  v'alues  of  the 
decay  distance  D  are  given  in  nautical  miles.  The  short  period  waves 
disappear  rapidly,  leaving  only  the  longer  period  waves  to  decay  mere 
slowly. 


The  areas  under  each  spectrum  are  related  to  the  square  of  the 
corresponding  significant  wave  height,  and  the  modal  period  is  clc  -e 
to  the  significant  wave  period.  Figure  34  shows  the  change  in  sigr  .fitani 
wave  height  and  period  as  a  function  of  decay  distance.  P’rom  this  iigure 
it  is  seen  that  most  of  the  decay  takes  place  during  the  first  1000  miles 
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FIGURE  32  FORECASTING  CURVES  FOR  WAVE  DECAY 
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FIGURE  33  DECAY  OF  WAVE 

SPECTRA  WITH  DISTANCE 
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after  which  little  change  occurs,  which  is  in  agreement  with  the  idea 
proposed  by  Munk(1964).  The  relationships  of  Figure  34,  having  been 
ba^ed  on  empirical  data,  do  not  indicate  how  the  energy  is  dispersed, 
dissipated  or  transformed  from  one  frequency  to  another. 

No  relationships  have  been  developed  tc»  -ht.w  how  the  wave 
spectrum  decays  with  respect  to  time  at  a  particular  location.  The 
wave  spectra  data  of  Ijima  (1957)  show  that  the  shift  in  modal  pericd  is 
opposite  to  that  given  in  Figure  33,  as  should  be  expected.  Figure  55 
represents  frequency  spectrum  for  increase  cl  waves  and  for  dec^ea-t- 
of  wave.s,  reprodvced  from  Ijima  (1957).  The  growth  rf  wave  energy  is 
from  high  frequency  (lew  period)  and  the  decay  cl  wsve  energy  is  f-.  m 
lew  frequency  (Icmg  period). 

The  c'.mbinaticn  of  wind  waves  and  swei.  i  '-.r-  lead  tc  a 
peaked  spectrum.  This  is  particularly  true  f:T  h  jrricare  conditionr 
where  the  hurricane  wiiids  act  at  nearly  a  right  ^r.gle  to  the  direction 
of  swell  propagated  ahead  of  the  storm  as  the  storm  moves  cover  the  sweK. 
A  typical  illustration  of  these  conditions  is  shown  for  Hurricane  Donna 
by  Bret  Schneider  (1961)  in  an  appendix  to  the  paper  by  Caldwell  and 
Williams  ( 1961)  as  published  in  Ocean  Wave  Spectra  ( 1965).  A  demcn- 
stratinr  of  this  phenomenon  is  shown  illustratively  in  Figure  36  for  the 
period  spectrum  and  which,  when  transformed  to  the  frequency  spectrum, 
is  shown  in  Figu'e  37,  remembering  that  S(f)  S(T).  It  might  be 

noticed  that  a  careful  inspection  of  Figures  56  and  3  7  reveals  the  fac  t  that 
for  certain  applications  the  period  spectrum  ha-  certain  advantages  ever 
the  frequency  spectrum.  It  appears  that  if  the  frequency  spectrum  were 
plotted  on  a  log  scale,  there  would  be  a  clarification  of  the  apparer.* 
difflc  ult  ies . 

Figure  32  Is  based  c*ntirely  on  meager  rmp’:-.’coi  data.  Tbf  -- 
"claticn -^hips  take  into  account  the  generating  ’t'.cb  -nci  s  i 

TesK-t  give  a  family  of  relationships  for  wave  de>'^>,  wheifc-is  •  r  ^ 
work  '  f  Sverdrup  and  Marti’,  (1947)  gives  a  single-  r  el  r*  ’ ■:  r.  -  hip  ’’  w  ‘'  • 
dvcA  v.  It  should  be  mentioned  that  the  wc  rk  -f  Pie  r--  r.-,  iSe  amar  r  ^tc 
3ame-  (  1955),  as  outlined  in  Hydrographic  Office  Public  ati-3n  b')  ' ,  c-itex - 
•ihfc  he  =t  scientific  means  for  describing  the  decay  waves.  Farther 
research  and  collection  of  field  data  is  necessary  befere  the  lij^e.'’ 
te  'hniques  proposed  in  H.  O,  publication  603  can  be  inrpiemented  i-'t 
I  Meal  applications. 
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WAVE  ENERGY  SPECTRUM  IN  THE 
BUILD-UP  AND  DECAY  OF  WAVES 
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FIGURE  36  EXAMPLE  OF  PERIOD 

SPECTRA  OF  COMBINED 
LOCAL  STORM  AND  SWELL 
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FIGURE  37  EXAMPLE  OF  FREQUENCY  SPECTRUM  OF  COMBINED  LOCAL 

STORM  AND  SWELL 


6500 


VI. 


WAVE  STATISTICS 


Wave  statistics  are  defined  in  terms  of  probability  of  occurrence 
or  recurrence  intervals,  for  example  the  average  number  of  years 
required  for  a  particular  value  of  the  significant  wave  height  to  be 
equalled  or  exceeded.  This  definition  does  not  necessarily  have  the  same 
statistical  meaning  as  wave  variability.  Wave  variability  is  reserved  for 
wave  height  or  period  distribution  for  a  particular  continuous  wave  record. 

In  this  section  probability  is  defined  as  the  percent  of  time  a 
particular  event  is  expected  to  occur.  The  cumulative  probability  is 
the  percent  of  time  a  particular  event  is  equalled  or  exceeded.  The 
recurrence  interval  is  the  time  required  for  a  particular  event  to  recur. 

In  order  to  determine  cumulative  probability  and  recurrence  intervals 
of  wave  heights  from  severe  storms  which  might  be  prognosticated  from 
climatological  events,  several  methods  of  approach  might  be  utilized. 

If  data  were  measured  over  a  sufficient  length  of  time,  this  would  be 
useful.  Otherwise,  wave  hindcasts  from  past  meteorological  weather 
maps  can  be  made  for  presentation  of  the  data.  The  method  of  Beard 
(19'j2)  has  been  used  successfully  for  small  samples  --  20  to  40  years 
of  records  »-  by  hydrologists  and  hydraulic  engineers  to  predict  recur¬ 
rence  intervals  for  peak  floods  in  engineering  studies  of  watersheds, 
reservoir  capacity  and  dam  construction. 

The  equations  from  Beard  (1952)  for  determining  probabilities 
and  recurrence  intervals  are  as  follows: 


and 


P  = 


100 


lOOY 


(86) 


where  S 


total  number  of  occurrences  on  record. 


s  =  the  summation  of  occurrences,  beginning  with  the  lowest 
value  to  any  successive  higher  value  until  s  =  S. 

P  -  the  cumulative  probability  that  an  event  is  equal  to  or 
less  than  a  particular  value. 

Y  -  the  number  of  years  of  record. 

I  -  the  recurrence  interval  in  years  corresponding  to  the 
probability  P  and  the  number  of  years  of  record  Y  . 
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Generally,  the  probability  P  versus  the  magnitude  of  the  v-veM 
is  aplotted  on  probability  paper  and  a  smooth  curve  is  constructed.  Fr  m 
the  smooth  curve  values  of  the  event  and  the  corresponding  prcbab'L/y 
are  determined,  and  the  recurrence  interval  I  is  then  determined  Irt-m 
Eq.  (87). 


Instead  of  using  the  graphical  plot  on  probability  paper,  c  ne  couiO 
use  to  a  better  degree  of  accuracy  the  Weibull  distribution  funrticr,  Eq. 
(33): 


Ji  n 


Y 


:n  [r[i  +  _L) 


:sH' 


If  the  data  fallow  a  nearly  linear  relationship  according  tc  the 
above  equations,  then  a  graphic^  determination  can  be  made.  If 
required,  a  statistical  least  squares  technique  can  be  used  to  obtain  the 
best  fit  parameters  A  and  M  .  Once  the  proper  analytical  equation 
has  been  determined  for  the  probability  P  ,  then  the  recurrence  interv-ii 
I  can  be  determined  as  before. 


Jasper  (1956)  presented  statistical  distributions  of  significar.t 
wave  heights  based  on  about  6-1/2  or  6  years  of  Weather  Burea’;  data 
for  certain  stations  of  the  North  Atlantic.  These  data  resulted  in  a 
linear  relationship  on  log  normal  probability  paper.  Figure  ;s8  ir 
reproduced  from  Jasper  (1956). 

More  recent  data  on  wave  statistics  have  been  pre.-er.red  :r  'ab  ■ 
form  by  M.  Darbyshire  (1963).  These  data  are  in  term>  ci  *he  movimurr 
wave  measured  and  the  corresponding  wave  period. 
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PROBABILITY  OF  EXCEEDING  WAVE  HEIGHT  ( PERCENT  ) 
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FIGURE  38  COMPARISON  OF 
WAVE  HEIGHT  DISTRIBUTIONS  DERIVED 
FROM  VISUAL  OBSERVATIONS  AND  FROM 
MEASUREMENTS  OF  WAVE  HEIGHTS  AT 
ATLANTIC  OCEAN  STATIONS  I  AND  J 
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PROBABILITY  OF  NOT  EXCEEDING  WAVE  HEIGHT  (PERCENT) 


vn.  WIND  SPEED  VERSUS  WIND  SPEED 


This  section  on  wind  speed  versus  wind  speed  is  indeed  a  subject 
about  which  much  discussion  can  be  promoted.  Although  there  are  con> 
troversies  about  use  of  the  proper  wind  speed  in  wave  forecasting,  there 
can  be  no  disagreement  that  much  research  is  required  to  standardize 
wind  speed  as  related  to  wave  generation.  It  appears  that  each  author 
of  a  particular  wave  forecasting  technique  has  his  own  definition  of  wind 
speed,  whereas  the  user  of  the  wave  forecasting  techniques  sometimes 
uses  a  mixed  definition.  It  is  no  wonder  that  the  various  methods  have 
been  ciritcized  by  the  users. 

Pierson  (1963)  for  the  first  time  gave  a  reminder  that  the  wind 
data  obtained  by  Neumann  (1948  and  1952)  was  actually  measured  at  7.  5 
meters  above  the  sea,  aboard  a  ship,  which  itself  was  bobbing  in  the  sea. 

A  hand  anemometer  extended  on  a  pole  away  from  the  ship  was  used  to 
record  the  wind  speed,  and,  at  the  same  time,  the  waves  were  observed. 
The  type  of  average  wind  speed  determined  is  difficult  to  know.  Perhaps 
some  of  the  scatter  of  the  data  given  in  Figure  18,  reproduced  from 
Neumann  (1952),  is  due  to  conditions  under  which  the  measurements  were 
made.  However,  the  averages  can  still  be  of  importance. 

An  attempt  by  Pierson  (1963)  was  made  to  interpret  wave  spectra 
in  terms  of  the  wind  profile  instead  of  wind  measured  at  a  constant 
elevation,  and  Pierson  and  Moskowitz  (1963)  proposed  a  wave  spectra 
based  on  wind  data  measured  at  a  constant  height  of  19.  5  meters  above 
mean  sea  level.  This  was  done  to  avoid  using  drag  coefficients  by  other 
authors  for  reducing  the  wind  speed  to  some  other  standard  elevation. 
Pierson  (1963)  also  attempted  to  reconcile  the  differences  between  the 
various  spectra  on  the  basis  of  the  elevation  at  which  winds  were  me&.sured, 
and  he  advocated  longer  averages  for  wind  speed  in  order  to  obtain  better 
correlation  with  the  averaged  wave  data. 

Darbyshire  (1952,  1955)  related  wave  data  with  the  geostrophic 
wind  speed  calculated  from  the  surface  pressure  gradient.  Since  the 
isobars  were  smoothed  through  the  recorded  pressure  data,  the  geo¬ 
strophic  wind  represented  some  type  of  a  longer  time  average. 

Bretschneider  (1951,  1959)  calculated  a  surface  wind  speed  from 
the  geostrophic  wind  speed,  taking  into  account  curvature  and  air -sea 
temperature  differences.  The  corrections  due  to  curvature  and  air-sea 
temperature  differences  were  based  on  the  original  empirical  Scripps 
data  given  by  Arthur  (1947).  These  data  were  correlated  with  the  ob¬ 
served  or  recorded  one -minute  average  wind  speeds.  Figure  39, 
based  on  the  geostrophic  wind  duration,  and  figure  40,  based  on  the 
Scripps  data  of  Arthur  (1947),  can  be  used  to  calculate  a  mean  surface 
wind  speed,  presumably  for  an  elevation  of  10  meters  above  the  mean 
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5mb  ISOBAR  SPACING  -  DEGREES  LATITUDE 
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FIGURE  39  GEOSTROPHIC  WIND  SCALE 
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3mn  ISOBAR  SPACING  -  DEGREES  LATITUDE 


RATIO  OF  SURFACE  WIND  TO  GEOSTROPHIC  WIND  (U/V^) 
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FIGURE  40  SURFACE  WIND  SCALE 
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sea  surface.  Sometimes  the  pressure  gradients  are  not  too  well  defined, 
particularly  for  lighter  winds,  but  in  any  case  the  calculated  wind  speeds 
should  agree  with  the  measured  wind  speeds  to  within  t  15%. 

The  data  of  Goodyear  (1963)  show  a  considerable  amount  of 
scatter  between  the  calculated  wind  speed  and  the  observed  wind  speed, 
and  states  that  the  "calculated  wind  speeds  represent  a  10  to  15  minute 
average"  whereas  the  measured  wind  speeds  are  at  the  most  several 
minutes  duration.  Figure  41,  reproduced  from  Goodyear  (1963),  shows 
the  deviations  between  observed  one-  or  two-minute  average  wind  speed 
and  the  calculated  10  to  15  minute  average  wind  speed  from  pressure 
gradients.  It  is  no  wonder  that  Wilson  (1963)  could  not  calculate  wind 
speeds  from  the  pressure  gradients  and  obtain  a  correlation  with  the 
reported  wind  speeds.  Instead  Wilson  (1963)  uses  the  reported  wind  speeds 
to  calculate  isolines  of  constant  wind  speed,  but  this  must  also  entail  an 
averaging  process,  an  average  with  respect  to  distance.  Further,  when 
using  the  space -time  wind  field  concept  there  is  still  another  average: 
a  time  average  from  one  weather  map  to  the  next.  Some  of  the  ship 
reports  of  winds  are  winds  measured  at  elevations  different  from  19.  5 
meters  above  the  sea.  One  can  hardly  help  wonder  what  wind  speed 
is  being  talked  about. 

Certainly  there  are  directional  spectra  of  winds  (which  vary  with 
elevation),  just  as  there  are  directional  spectra  of  waves.  Horizontal 
variation  in  wind  speed  and  direction  must  be  considered  of  great  im¬ 
portance  in  wave  generation,  particularly  when  one  discusses  effective 
fetches  and  durations  as  well  as  fully -developed  seas.  Until  more  care 
is  taken  in  reporting  winds  and  defining  wind  variability,  one  cannot 
expect  any  novel  advances  in  wave  forecasting.  Even  then  it  will  be 
required  to  correlate  new  wave  spectra  data  with  new  wind  speed  deter¬ 
minations. 

Presently  the  forecaster  must  use  the  existing  tools  for  calculating 
waves.  He  must  use  the  type  of  wind  analysis  for  making  the  forecast 
which  is  compatible  to  the  methods  used  to  develop  the  forecasting 
relationships  and  procedures.  A  change  in  the  procedures  on  which  the 
tools  are  originally  based  should  also  require  a  change  in  the  tools,  or 
else  the  accuracy  of  the  forecasts  will  become  worse  instead  of  better. 
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OBSERVED  WIND  SPEED  (M«) 


COMPUTED  WIND  SPEED  (kts) 

(  AFTER  GOODYEAR ,  1963  ) 


FIGURE  41  COMPUTED  vs.  OBSERVED 

SURFACE  WIND  SPEED  FOR  43 
RANDOMLY  SELECTED  POINTS 
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